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1. Introduction 



An interesting solution to the equations of motion of two dimensional string theory is 
the semi-infinite cigar Jl|J^,|[, described by the metric and dilaton 

ds 2 = k \dr 2 + tanh 2 r d9 2 ] 

(1-1) 

$ — $o = — 2 log cosh r 

where 9 is a periodic coordinate, # ~ # + 2tt, labeling the location around the cigar, and 
r > is the direction along the cigar, with r = corresponding to the tip (see fig. 1). The 
string coupling depends on r - it goes to zero far from the tip of the cigar, and attains 
its maximal value, g s = e*°, at the tip. k is a free parameter which governs the overall 
size of the cigar. 
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Fig.l. The semi-infinite cigar. 



Some of the contexts in which the geometry ( |1 . 1| ) appears in string theory are: 
^1) This geometry can be thought of as a Euclidean version of a 1+1 dimensional black hole 



||. The Minkowski black hole can be obtained from (|1 . 1| ) by analytically continuing 
6 — > it. Furthermore, string propagation in the geometry ( |1 . 1| ) can be described 
by a coset conformal field theory; the Minkowski version corresponds to the coset 
SL(2, R)/U(l), while the Euclidean one is su(2) /^(-O - The parameter k in ( |1 . 1[ ) is 
the level of the SL(2) current algebra that enters the construction;!! the mass of the 
black hole, M, is encoded in the value of the string coupling at the tip of the cigar, 
M oc exp(-2$ )- 

(2) The near-horizon geometry of non-extremal iVS'S-branes is a tensor product of the two 
dimensional black hole geometry (|1 . 1| ) , R 5 (the spatial directions along the fivebranes) 
and S 3 (the angular directions transverse to the fivebranes) ||. The corresponding 
worldsheet theory is a supersymmetric version of the above coset CFT. In this case k 
is the number of iVS'S-branes, and the string coupling at the tip of the cigar measures 
the energy density of the fivebranes (see M for a recent discussion). 



1 The background ( fL.lD is in fact only valid for large k. For finite k there are corrections that 
were found in [01. 
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(3) String theory on a Calabi-Yau manifold near a point in moduli space at which the 
manifold develops an isolated singularity is described by a geometry that includes ( |TTT| ) 
as a factor f7||| . The level k depends in this case on the particular singularity, and the 
string coupling at the tip of the cigar corresponds to a non-perturbative energy scale 
in the theory - the energy of D-branes wrapped around cycles whose size vanishes at 
the singularity. 



Further study of the geometry (|1.1|) would thus be useful for understanding the physics 
associated with horizons and singularities in string theory. In particular, one would like 
to obtain a better understanding of the thermodynamics and statistical mechanics of the 
corresponding black holes (or black branes). 

A promising framework in this regard is holography. In all three examples mentioned 
above, string theory in the geometry ( |1 . 1| ) is believed to be holographically dual to a non- 
gravitational theory. In the case (1) the dual is the "old matrix model" - Matrix Quantum 
Mechanics (MQM) in the double scaling limit - while in the other two examples it is Little 
String Theory (LST). 

In this paper we will focus on the first example, the 1 + 1 dimensional black hole. We 
will propose a solvable formulation of two dimensional string theory in the background 
( |1 . 1| ) , based on MQM. This approach to two dimensional string theory in flat spacetime 
with a condensate of the tachyon field |J led to a solution of the theory to all orders 



in string perturbation theory ||T0| , |rT| , |12l . |T3jl ; see |]l^ , |15| , |T6|| for reviews. The collective field 
theory of the singlet sector of MQM was identified with the dynamics of the tachyon field 
in this background |fL7l| . Some early attempts to find a matrix realization (in the singlet 
sector of the matrix model) of two dimensional string theory in the background ( |1 . 1| ) are 
described in 0,0^3 . 

As mentioned above, one of the main motivations for looking for a holographic de- 
scription is to study the statistical mechanics of these black holes. Standard black hole 
thermodynamics arguments imply that there is a Hagedorn density of black holes: 

p{M) ~ M Sl e^ M . (1.2) 

The inverse temperature 0h can be read off the geometry ( p. . 1|) , while s\ is given by a one 
loop calculation in string theory. The prediction ( |1.2|) is particularly interesting since, at 
least perturbatively, two dimensional string theory has very few states (one field theoretic 
degree of freedom in 1 + 1 dimensions). Thus, it is interesting to verify the thermodynamic 
prediction ( |1.2| ) by a direct counting of states. In other examples, this proved easier to do 
in the holographically dual picture, and one might hope that the same will happen here. 

The starting point of our analysis is an observation made by V. Fateev, A. Zamolod- 
chikov and Al. Zamolodchikov [^T|. These authors conjectured that the SL(2)/U(1) coset 



CFT is equivalent to the Sine-Liouville model, i.e. c = 1 CFT coupled to a Liouville field, 
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with the cosmo logical constant tuned to zero and the scale set by a winding mode of the 
c = 1 field. This correspondence has not been proven but evidence for its validity has 
been presentedi |2l[] . Our main purpose here will be to define and study the Sine-Liouville 
background of two dimensional string theory using 

From the point of view of the usual Liouville background of two dimensional string 
theory, the Sine-Liouville vacuum corresponds to condensation of vortices (winding modes) . 
We will treat the vortices as a perturbation, and will show that in the cases of interest, 
the process of taking the strength of the perturbation to infinity (corresponding to vortex 
condensation) can be studied using MQM. Some steps in this direction were taken in the 
past. In particular, G. Moore |2t| studied the Sine-Gordon model coupled to gravity by 



treating the Sine-Gordon term as a perturbation. The resulting structure was further 



discussed in []24j where it was shown that the model describes a renormalization group 
flow from the c = 1 model coupled to gravity in the UV, to a set of decoupled c = models 
coupled to gravity in the IR. As we will see below, in a certain region of the parameter 
space of the model, there is another possible critical point - the Sine-Liouville model. 

In the matrix model formulation, it was conjectured [p5| , |26[| and then demonstrated 



||27|| that vortices on the worldsheet correspond to U (N) non-singlet states in MQM. Their 
wave functions are described by Young tableaux of U(N); the vortex charge is equal 
to the number of boxes. It is also possible to study vortices, at least in the spherical 
approximation, using the dual matrix description based on discrete time HTT| , |26| , p8| . 

By using the connection between vortices and non-singlet states in MQM, we will 
construct an integrable system, the infinite Toda chain hierarchy, that interpolates between 
the usual c = 1 string and Sine-Liouville (or two dimensional black hole ||21[| ) backgrounds 
structure allows one to compute the partition sum and correlation functions of the theory 
to all orders in string perturbation theory. As an example we will give explicit expressions 
for the partition functions on the sphere and on the torus@. 

The plan of the paper is as follows. In section 2 we review the precise statement 
of the correspondence between the SL(2)/U(1) coset CFT and the Sine-Liouville model, 



and some of the evidence for it. In section 3 we use the results of |23|,^3| to study the 
c = 1 string theory (at tree level) in the presence of the Sine-Liouville interaction. We 
show that in a certain range of the parameters defining the model, there is an obstruction 



2 The Sine-Liouville model is the massless limit of the "sausage model" considered in 

3 A direct argument in favour of the description of the 2d black hole in terms of compact 
MQM, based on the similar conjecture relating thermal YM to the ADS black hole, was given to 
us by A. Polyakov in 1998. 

4 Toda integrable structure of the c = 1 string theory perturbed by purely tachyon source has 
been first discovered by Dijgraaf, Moore and Plesser in |29[ |. Our description is related to their by 
T-duality. 
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to turning off the cosmoiogical constant and using the Sine-Liouville interaction to set 
the scale (associated with an instability of Sine-Liouville theory), but in the range of 
parameters relevant for the correspondence of section 2 this obstruction disappears and 
the Sine-Liouville/black hole phase is stable. 

In section 4 we introduce the MQM describing vortices and reduce it to a one-matrix 
model whose degrees of freedom correspond to the matrix holonomy factor around the 
compactification circle. We find that the partition function of this model is a r-function 
of the infinite Toda chain hierarchy! We prove this statement for a more general partition 
function, which can be thought of as the generating function of all (multi-) vortex correlators 
in two dimensional string theory. 

In Section 5 we use the Toda chain equations and the KPZ-DDK scaling properties 
of the partition function to obtain ordinary differential equations for the genus h partition 
functions. We solve these equations explicitly for the sphere (h = 0), thus proving a 
conjecture of [f23|, and for the torus (h = 1) which is a new result. 

In section 6 we discuss the implications of our results for black hole thermodynamics. 
In particular we compute the coefficient a in ( |1.2| ) and show that the Hagedorn temperature 
is associated with a phase transition. In section 7 we briefly describe the thermodynamics 
of black holes in fermionic two dimensional string theory, which helps to clarify some 
aspects of the bosonic analysis. In section 8 we comment on our results and describe some 
open problems. Some of the technical details are contained in the appendices. 



2. On the equivalence of the cigar and Sine-Liouville CFT's 

As mentioned above, the Euclidean black hole worldsheet CFT ( |1 . 1| ) corresponds to 
the coset 

H t ,U W ; H t S (2.1) 

It can be studied using standard coset CFT techniques, and is closely related to CFT on 
Euclidean AdSs(= H^); see e.g. |30[| . Below we list some properties of this theory. 

The central charge of (|2.1| ) is given in terms of the level of the underlying SL{2) 
current algebra, k, by 

3k , . 

1. (2.2) 



k-2 

In particular, for k = 9/4, c = 26, and the coset is a good classical solution of two dimen- 
sional bosonic string theory. An important set of observables corresponds to momentum 
and winding modes on the cigar, Vj- m> fn- These operators have scaling dimensions 

a j(j + l) m 2 - + 1) f to 2 

^j;m,m — fe — 2 ' ■ 7 ' ;m ' m ~~ k - 2 k ' 



The fact that the generating function for the tachyon amplitudes in the c = 1 string theory 



posesses Toda lattice symmetry has been established in [29]. 
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where m, m are given by 

m=^(ni+n 2 fc); rh = -^(n 1 - n 2 k); m,n 2 e Z. (2.4) 

One can think of the integers n\ and ni as momentum and winding around the cigar (in 
the 9 direction in figure 1). As is clear from the geometry of the cigar, the momentum, 
m — m, is conserved in this CFT, while the winding, m + m, is not. 

Spherical two and three point functions of the observables Vj ;Tri) fn were computed in 
|2T],|3D[]. The two point function is (suppressing the dependence on the separation of the 
worldsheet locations of the vertex operators) 

„ . , . r(i - f±i)r(-2j - i)r(j -m + i)r(i + j + m) 



r(f^)r(2j + 2)T(-j - m)T{m - j) 

(2.5) 

where 

^(l-j^) 



The result for the three point function is more involved and can be found in |3(J . 

String perturbation theory is an expansion in exp(— 2$). On the cigar, the dilaton $ 
varies as in (|1 . 1| ) , but string perturbation theory still makes sense - it can be thought of 
as an expansion in exp(2$o) 5 the value of the dilaton at the tip (which is proportional to 
1/M). Thus, the genus expansion on the cigar is a 1/M expansion. 



The dual Sine-Liouville theory |TlJ is described by the Lagrangian (we set a' = 1) 
1 



L 

A-k 



{dxf + (d(f>) 2 + QR(/) + \e H cosR{x L - x R ) . (2.7) 



The central charge of this theory is c = 1 + 6Q + 1; comparing to (|2.2|) we find that 



« 2 = ih- (2 - 8) 

The radius of x is R = y/k, the same as the radius of 9 far from the tip of the cigar in 
( pL.l|) . The cosine interaction in (|2.7|) is thus the lowest lying winding mode, with winding 
number equal to one. Requiring that the scaling dimension of the sine-Liouville interaction 
is equal to one, leads to 

\R 2 -\b(b + 2Q) = l (2.9) 

with the solution 

1 



- = -y/k=2. (2.10) 
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Very far from the tip of the cigar, (LI) describes a cylinder with a dilaton field that varies 
linearly along it. The same is true for the geometry corresponding to (|2.7| ) far from the 
potential wall (i.e. at <fi — > oo). The two cylinders can thus be identified (for large k) via 



a % (2.11) 

v ~ — — . 

Vk 



The non-trivial statement of |2T[] is that the two CFT's also agree when one includes the 
interaction that cuts off the large string coupling region (j> — > — oo. This is far from obvious, 
since in one case the strong coupling region is eliminated by changing the topology of the 
cylinder to that of the cigar (|1 . 1| ) , while in the other this is achieved by turning on the 
potential in (|2.7| ). 

The relation between the coset and Sine-Liouville CFT's is a strong-weak coupling 
duality on the worldhseet. The cigar CFT becomes weakly coupled in the limit k — > oo, 
where the geometry (|1 . 1|) becomes weakly curved. In that region the wavefunction of the 
Sine-Liouville potential, which has the large <fi behavior 

#(0)~e (Q -^ (2.12) 

goes rapidly to zero as — > oo, i.e. ^ is supported in the region of small <p and the theory 
is strongly coupled. The semiclassical limit in Sine-Liouville is Q — > oo (i.e. k — > 2 ( |2.8| )). 
In this case the wavefunction ( |2.12| ) is supported at (p — > oo, far from the potential wall 
( |2.7| ) and the theory is weakly coupled. The cigar is highly curved in this limit; thus the 
coset CFT is strongly coupled. 

The operators Vj- mj rh mentioned above are described in the Sine-Liouville model a 

Vj.^rn <- jP^+iWa+M, (2.13) 

where 

PL — —=r + n 2 R 

ti 

PR = ^-n 2 R (2.14) 
P = 2Qj. 

Indeed, the scaling dimensions 

A = \A - \P{P + 2Q); A = \p 2 R - \p(P + 2Q) (2.15) 



More precisely, this is the form of the vertex operators as (f> — > oo. 
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agree with ( |2.3|) . Eq. ( |2.14 ) also makes it clear that ni and n<i are the momentum and 



winding around the cylinder (or cigar). Note also that like in the cigar CFT, in (|2.7|) the 



momentum p^ + Pr is conserved, while the winding number p^ — Pr is broken. In the 
cigar CFT the reason for that is that winding can slip off the tip of the cigar; in Q2.7| ) the 
interaction breaks this symmetry explicitly. 

The Sine-Liouville model (|2.7| ) has similar scaling properties to Liouville theory. For 
example, the standard scaling analysis implies that the partition sum has the following 
perturbative expansion: 

T{\g s ) = Y,^h{g s \~^) , (2.16) 

7i=0 

where Th is the genus h partition sum, and g s is the string coupling constant.0 As in 
Liouville theory, the physics depends only on the combination ^A" 7 ^, which can be 



thought of as the effective string coupling in the background (|2.7| ). From now on, we will 
absorb g s into A (i.e. set g s = 1). 

We see that the string coupling expansion in Sine-Liouville theory is a large A expan- 
sion, like in ordinary Liouville theory. In string theory on the cigar, as mentioned above, 
the string coupling at the tip of the cigar is related to the mass of the black hole via 
gl ~ 1/M where M is measured in Planck units. Therefore, in the correspondence of [21 
A is related to the mass of the black hole, M, 



2 



M^\—. (2.17) 

An important piece of evidence for the equivalence of the cigar and Sine-Liouville CFT's 
is the agreement of a large class of two and three point functions on the sphere in the two 



models [|21] . We will next review this agreement for the two point function; the three point 
functions can be discussed along the same lines. 

The two point functions of Vj ;m> fh in the coset CFT, given by (|2.5| ), exhibit a series 
of poles. As explained in || (for the super symmetric similar analysis applies to 

the bosonic case discussed here), the poles of the first V function in the numerator play 
a different role from those of the rest of the factors. These poles occur at values of j for 
which 

2 ?' + 1 , 
l-J—^ = -n; n = 0,l,2,---. (2.18) 

Near such a pole, the amplitude has the form 



7 defined as the value of exp <E> at some arbitrary point along the cylinder. 
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(Y, 



I : in,mVj;—m, — m) ~ 



(k - 2)[v(h)] 2 3 +1 (-l) n T(-2j - 1)T(J -m+ 1)T(1 + j + fh) 



n!(n+l-f±± 



) r(fS) r ( 2 i + 2 ) r (-^' - ™) r (m - 3) 



(2.19) 

As observed in |2"T{] , the poles ( ^.18| ) are also obtained in Sine-Liouville theory (|2.7|), and 
the residues agree with (|2.19| ). We next review this calculation. 
The two point function (|2.5|) corresponds in Sine-Liouville to 



A{pl,Pr) = (e 



iPLXL+iPRXR+P(l> e -iPLXL-iPRXR+l3<j>\ 



(2.20) 



(see ( |2.13| ), ( |2.14p ). The theory ( |2.7| ) is interacting, and in general such correlation func- 
tions are difficult to compute (although much progress has been made on this problem 
in the last ten years). However, in some cases one can compute the residues of certain 
poles (which will turn out to correspond precisely to ( |2.18| )) as follows (see e.g. []32| for a 
review). Split the path integral over the Liouville field <p into a zero mode integral and a 
path integral over the non-zero modes. Performing the integral over the zero mode, one 
formally finds (after absorbing an overall constant into the definition of the path integral) 



ip L x L +ip R x R +l34> e -ip L x L 



-ip R XR+P4> 



e bc ^ cos R(xl — xr) 



)a=o- 



A(p L , PR ) = \ s T(-s)(e* 

U J (2-21) 

The expectation value in (|2.21|) is understood to exclude the zero mode of (f>; note also 
that it is performed with the free action (A = in ( |2.7|) ). s is the KPZ scaling parameter, 
which using ( |2~TU| ), (|2~TJ) is 

2(2j + 1) 



k-2 



(2.22) 



In general the representation ( |2.21 ) is highly formal since it is difficult to make sense 
of the non-integer power of the interaction in the correlator. In situations where s is a 
non-negative integer, the correlator in ( |2.21| ) does seem to make sense but the prefactor 
diverges. The nature of the divergence is well understood [32j. Amplitudes with s G Z + 
are bulk amplitudes - they correspond to processes that can occur anywhere in the infinite 
region far from the Sine-Liouville wall and the divergence in question is nothing but the 
volume of that region.! Regularizing the theory by cutting off the region <fi — » oo (far from 
the wall ( |2.7[ )) amounts to replacing 



X s T(-s) 



'-l) s+1 A 
1 s\ A' 



(2.23) 



Note that these amplitudes are in general sensitive to the structure of the wall, since they in- 
volve bulk interactions between the incoming quanta and those that make up the wall. Amplitudes 
with s = are insensitive to the structure of the wall, a fact that will play a role later. 
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where A is a UV cutoff. 

To summarize, the residues of the poies of r(— s) in ( |2.21| ) are given by bulk am- 
plitudes, which have an integral representation of Shapiro- Virasoro type. These can be 
computed by standard techniques. 

When s is odd, the bulk amplitude vanishes. The reason is that winding number, 
which is broken by the Sine-Liouville interaction, is actually a good symmetry in the bulk 
of space, and for s G 2Z + 1 the correlator ( |2.21| ) does not conserve winding. Thus, only 
poles with 

s = 2(n + l); n = 0,1,2,- ■■ (2.24) 

survive (the amplitude with s = is special and can be shown to vanish in this case). 
Comparing to ( |2.18| ) we see that the amplitude A(pl,Pr) (|2.21[) has poles at the same 
places as (|2.5|) (the rest of the poles of (|2.5|) do not correspond to bulk physics; they have 
to do with the existence of certain normalizable states living near the tip of the cigar). 
Near the poles, ( j2.21| ) is given by 

A 2(n+1) 



2 n+i( 2n + 2)!(n + l)l 2 (n + 1-M±± 



l e ipLX L +iPRXR+l3(t> e -ipLXL-'iPRXR+P4> 



JR(x L -x R ) 



n+1 



-iR(x L -x R ) 



n+1 



(2.25) 

As usual, all but three of the vertex operators in (|2.25|) are integrated over the complex 
plane. These integrals can be evaluated (they are generalizations of those discussed e.g. in 
||32j|), and shown to agree with ( |2.19|) . We will not discuss the details of these calculations 
here. 

One might worry that this agreement is coincidental since one can rescale the op- 
erators Vj-m^fh by a (j;m, m) dependent factor. However, the relative normalization of 



the operators entering the correspondence ( p. 13 ) is in fact fixed by requiring that their 
wave- functions coincide in the region very far from the tip of the cigar or Sine-Liouville 
potential {e.g. they can be taken to be incoming waves of unit strength). This makes the 
comparison of the two-point functions, which upon analytic continuation to Minkowski 
spacetime correspond to the scattering matrices off the tip of the cigar or Sine-Liouville 
potential, meaningful. The authors of J2l| also showed that the three point functions of 
the operators (|2.13|) exhibit similar poles, and the residues again agree. 

It is useful to mention for future reference that the correspondence between the cigar 
and Sine-Liouville conformal field theories has a (worldsheet) supersymmetric generaliza- 
tion proposed in ||. It relates the A^ = 1 superconformal coset model SL{2) /U{1) to the 
N = 2 Liouville theory! (see e.g. for a description of A^ = 2 Liouville). 



9 It might seem surprising that an iV = 1 superconformal field theory can be equivalent to an 
N = 2 superconformal one. In fact, it turns out that the N = 1 supersymmetric SL(2)/U(1) coset 



has an accidental TV" 



2 supersymmetry. It is a special case of the Kazama-Suzuki construction 
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3. Sine-Liouville as a perturbation of the c = 1 string 

As discussed in section 2, one approach to studying the black hole background in two 
dimensional string theory is to study the dual Sine-Liouville theory fl2.7f) . This is what we 
will do below. 

In bosonic two dimensional string theory the central charge (|2.2j ) is c = 26, which 
determines the different parameters that appeared in section 2 to be /c = 9/4, Q = 2, 
R = 3/2. We would like to use the matrix model approach to analyze this model. An 
immediate problem that one runs into is that in the matrix model A is set to zero, and 
instead the Lagrangian contains the Liouville potential ^<pexp(— 20). One can attempt 
to define (|2.7|) by means of conformal perturbation theory, i.e. consider the theory with 
Lagrangian 

L = — \(dx) 2 + (d(/)) 2 + 2R(/) + ^e" 2 ^ + Ae^ -2 ^ cos R{x L - x R )} (3.1) 

and treat A as a perturbation. General arguments lead one to expect that the radius of 
convergence of the resulting series is finite. Since we are interested in a continuation to 
A — > oo (or equivalently \i — > 0), we are not guaranteed that the series will converge. 
Indeed, we will see that in some cases there is an obstruction to taking this limit - one 
encounters a singularity at finite (real) A. 

This is not unexpected; the classical Lagrangian ( |3.1|) with fx = has a potential 
which is not bounded from below. Thus, at least for small R, where the cosine is slowly 
varying and the semiclassical picture is reliable, we expect the limit \x —>■ not to exist. 
As R increases, quantum effects (on the worldsheet) become important, and the vacuum 
energy could in principle be renormalized in such a way that the limit exists. We will 
see below that this is indeed what happens. Fortunately, for R = 3/2, the value relevant 
for the correspondence with the cigar, there is no singularity and the Sine-Liouville model 
exists. Of course, this has to be the case if the conjecture described in section 2 is correct. 

Consider, for example, the partition sum of ( |3.1| ) on the sphere. Define 

j p — r ( p ) f e ip(x L -x R )+(\p\-2)<j>_ (3_ 2 ) 



r(i-p) 

After rescaling A appropriately, the partition sum corresponding to ( |3.1|) becomes 

^o(A ± ,^) = (e- A + r «- A - r -*>. (3.3) 

We are interested in the case A+ = A_ = -|A, but it will be convenient to keep A± 
independent. Expanding in A± and using winding number conservation, we find 

n=0 
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The coefficients in this series are 2n point functions on the sphere, which can in principle 
be computed in the matrix model. G. Moore [^3| conjectured a general form for these 
amplitudes based on an extrapolation of matrix model results from small n: 



(T2T\) = Rn\ fi 2 [(l - R)^- 2 ^ ^ 2 R \ % (n > 1) 
\ r -R/o p ii T(n(l-R) + 1) ~ ' 

The term with n = is the usual c = 1 partition sum on the sphere, 



(3.5) 



Fo(A = 0,^) = -^ 2 log^. (3.6) 
Combining it with ( ^.5| ) we find 

Fo(A, „) = -f f log f + £ 1 [d - fl)M«- 2 A + A_] " T ^-_%- + 2 l y (3.7) 



n=l 



Differentiating twice with respect to \i (and dropping a non-universal constant) we find 

r(n(2-fl)) 
(n(l-R) + l) 



Xo(A, M ) = ^Fo(A,^) = -iHog^ + it>]P ( ^ r ^^ ir (3-8) 



n=l 

where 



z = (-R- 1)A+A_^- 2 (3.9) 
measures the "dimensionless strength" of the perturbation A. The series in (^.8|) can be 



summed using equation 5.2.13.30 in 35 



^ Q n\T(n(a - 1) + b) 1 } 6-1 ' (1 - s) a ~ 1 ] 

Taking the limit 5 — * 1 and plugging in Q3.8Q (with a = 2 — R) leads to 

Xo(A,m) = -Rlogn + Rlog(l-s), (3.11) 

where s is given by 

For small A, expanding (|3.12|) one finds s ~ z + 0(z 2 ). The large z behavior depends on 



i?. Next we analyze this behavior function of R. 

We restrict attention to R < 2. The reason is that for R > 2 the A perturbation in 
( [3J]) qualitatively changes the behavior of the potential in the weak coupling (UV) region 
<fi — > oo. Instead of dying off, it now fluctuates with an amplitude that diverges as — > oo. 
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This is a reflection of the fact that in flat space (i.e. discarding </>), the A perturbation is 
non-renormalizable in this case. 

It turns out that there is also a difference between the behavior for R < 1 and R > 1. 
Consider first the case R < 1. The parameter z ( |3.9| ) is negative in this case and we would 
like to analyze the behavior as z — > — oo. However, looking at (|3.12|) we see that there is 
an obstruction: d s z is positive only for s > s c = —1/(1 — R). Thus, varying s along the 
negative real axis we can not probe z < z c = —(1 — R) 1 ~ R / (2 — R) 2 ~ R . 

The physics of this obstruction was explained in |[24)| . It corresponds to c = critical 
behavior which is expected to occur in the model (|3.1| ). The c = 1 field x (more precisely 
the T-dual field ) settles into one of the minima of the cosine potential and 

decouples, leaving behind a c = system (coupled to gravity). While this is more clearly 
visible in the fixed area representation (see |24|] for the details), it can also be seen directly 
from (|3.11|) , ( |3.12|) . Near the critical point we have z — z c ~ (s — s c ) 2 . Thus the leading 



singular behavior of ( |3.11|) as z — > z c is %o ~ ( z ~ z c)K which is the expected behavior for 
pure gravity. 

Thus it appears that for R < 1 the Sine-Liouville theory ( |2.7|) is unstable (as expected 
on general grounds - see the discussion after eq. ( |3.1| )) and can not be reached0 from 
(PHD by taking y. -> 0. 



We next turn to 1 < R < 2. In |24] it was shown that quantum corrections generate 
in this case a finite positive contribution to the effective cosmological constant. Therefore, 
one expects the limit \x — > to be non-singular. This is indeed visible in the formulae 
above. In this case z > ( |3.9|) and the relation between z and s ( |3.12j ) is monotonous. 
Thus, as z — > oo, s — ► 1. This is the critical point we are interested in, in which \i has 
been effectively switched off. As a check, note that in this limit the dependence of the 
correlator xo on A* drops out, as it should. Indeed, combining ( |3.11|) and ( |3.12|) and using 
we find the following algebraic equation for the susceptibility: 

ye^ + (R- l)A+A_e^ X0 = 1. (3.13) 

From here we find that at large A 



X0 ^^2 1Og[(jR_1)A+A - ] - (3 - 14) 



10 We note for completeness that we could discuss the system with positive z, in which case 
there is no obstruction to taking z — > oo. As is clear from ( |3.12 ), as z — > oo, s — > 1, and we seem 



to find a new fixed point. Positive z corresponds for R < 1 to imaginary A (see ( p.9[ )). Thus, for 
R < 1 one can study an analog of Sine-Liouville, but with an imaginary potential. This theory is 
non-unitary. Also, for R < 1 the winding number two perturbation is relevant and one can study 
flows involving it. 
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This is easily seen to be the correct KPZ scaling for the two point function of e _2< ^ in 
Sine-Liouville theory. 

One can in fact go further and find the whole small \i (or large A) expansion of the 
susceptibility.0 Namely, from ( |3.11| ) and ( |3.12| ) we obtain the following parametrization 
of the susceptibility 

Xo (A, y) = -R\ogy + R\ogt = -^-^ log ((R - 1)A+A_) + log(l - t), (3.15) 



where t(= 1 — s) is given by 



1 y (3.16) 



(1 - t)*=* 



and 



y = z-^ = ^ (3.17) 



((i2-l)A+A_)*=* 

Using the formula ( |3.10| ) again, we arrive at the following small \x expansion: 

R , tt v ^ x n , R (~y) n t K^r 



Xo (A, n) = %F (\, /*) = log ((R - 1)A + A_) + ^ £ ^ 



n + 1 



2-i? 

(3.18) 

The coefficient of // n /n! in ( |3.18| ) is the n point function of zero momentum tachyons in 
the Sine-Liouville background ( |2.7|) . According to the conjecture of plj] , for = 3/2 it 



should also be the n point function of these tachyons in the Euclidean black hole geometry. 

To summarize, we find that for 1 < R < 2, there is no obstruction to turning on a 
large (real) Sine-Liouville coupling (^.1|). There is a critical point corresponding to n = 0, 
A finite, which appears to be stable and according to the conjecture described in section 2 
is equivalent (for R = 3/2) to the black hole background Q2.1D . 



This raises an interesting question,0 what is the interpretation of the Sine-Liouville 
background with R ^ 3/2 in terms of the cigar geometry? Naively, changing the asymp- 
totic radius of the cigar creates a conical singularity at the tip, and thus breaks conformal 
invariance. One natural proposal is that changing the radius R in the Sine-Liouville de- 
scription corresponds in the cigar theory to changing the asymptotic radius in the same 
way (so the asymptotic geometries agree) and turning on a potential like that in (|3.1|). 
This might lead to a change of the radius of x as a function of <fi, such that the conical 
singularity at the tip is eliminated. It would be interesting to make this proposal more 
precise, but we will not attempt that here. 



li 



To do this it is useful to note that (3.13) is invariant under: 2 — R «-> y2 — R\ogfi 



^= \og[{R - 1)A+A_], ^xo - n^mXo- 
12 D.K. thanks S. Shenker and L. Susskind for a discussion of this issue. 
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In the next two sections we will construct a matrix model describing two dimen- 
sional string theory with the worldsheet Lagrangian (|3.1[) and use it to compute the string 



partition function. We will in fact study a more general problem corresponding to two 
dimensional string theory in an arbitrary winding background, 

L = ^-{{dxf + (d<j)) 2 + 2Rcf> + ^e' 2 * + y tne (HR-^^nR(x L - XR )^ (3 ig) 

The couplings t n carry winding number n; the model (|3.1|) , ( |3.3|) corresponds to the choice 

t n = 5 n ,iA + + <5 nj _iA_. (3.20) 

Of course, most of the couplings t n are in general irrelevant, which is reflected in the fact 
that the corresponding perturbations in (|3.19| ) are singular at <fi ~^ 00 • These couplings 
should be treated perturbatively, and the path integral ( |3.19|) should be thought of as 
the generating functional of correlation functions of these operators, as a function of the 
relevant and marginal couplings. 



4. A matrix model for the c = 1 string theory with vortices 

In this section we will describe a matrix model for the c = 1 string theory in the 
presence of vortices of arbitrary vorticity, or in continuum language the theory described 
by the Lagrangian ( |3.19| ) with t n 7^ 0. We will show that the grand-canonical partition 
function of this matrix model is a r-function of the Toda chain hierarchy and, as such, it 
satisfies the Toda lattice equation. This equation can be used to compute the partition 
sum of the model to all orders in the string coupling (or equivalently to all orders in 1/M 
in the black hole background). 

4-1. The role of twisted boundary conditions in MQM 

The usual c = 1 string corresponds to the fixed point in which all t n in (|3.19|) vanish 
and the scale is set by the cosmological constant \x. The (perturbative) dynamics is in this 
case well described by the singlet sector of the c = 1 matrix model and, as mentioned above, 
contains very few degrees of freedom. The Euclidean problem with finite R corresponds 
to finite temperature, T = l/2nR. When the temperature is sufficiently low (i? > 2 in 
our normalizations ( |2.14j ), ( |2.15| )), one can neglect the non-singlet states ||25|| . Above this 
critical temperature (R < 2) the vortices become important; studying the Lagrangian 
( |3.19| ) with t n 7^ is a useful way to incorporate them. 
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To describe the c=l string with vortex excitations (34.) we will use the twisted MQM 



introduced in | 27fl . Its canonical (fixed N) partition function is defined as a functional 
integral with respect to the one-dimensional N x N Hermitian matrix field Mj^{x) 

Z N (Q)= [ PM(x)e- tr r fid 4l^M) 2 +nM)] (41) 

JM(27ri?)=ntM(0)n 

where O is an arbitrary unitary matrix. Note that the path integral ( }4.1| ) depends on Q 
only through its eigenvalues, z\ = e ldl , • ■ ■ , zn = e l0N . The potential V can be taken to be 
any function with the appropriate critical behavior, e.g. 

V(M) = -M 2 g —M z . (4.2) 



The projection onto the singlet sector is achieved by imposing the twisted boundary con- 
ditions ( f4.1| ) and then integrating with respect to the twist matrix O. As we review 
in Appendix A, the singlet sector can be reduced to a system of free fermions corre- 
sponding to the eigenvalues (yi, ■ ■ ■ ,Un) of the matrix M, each compactified on the circle, 
yj(x + 2tvR) = yj{x) and moving in the potential V(y). The introduction of vortex exci- 
tations is achieved by including higher representations of SU(N) |25,27]. The eigenvalue 
dynamics in these representations is more complicated. 

To find the realization of the system ( |3.19| ) in MQM we next take a closer look at the 
role of Q in the large limit. The 1/N expansion of the free energy 

of the O-twisted matrix model ( f4.1| ) can be expressed in terms of connected planar graphs, 
embedded in the target-space circle. Namely 



T n (Q) = J2^n\^) N 2 ~ 2 \ (4.3) 



h>0 

where JF^^O) is the contribution of planar graphs with the topology of a sphere with h 
handles. We will use here the usual connection between planar graphs in MQM and string 
worldsheets: each planar graph E is viewed as a discretized worldsheet whose vertices, 
links and faces are labeled by v , £ and /, correspondingly, and the functional integral over 
the worldsheet fields (x, </>), where <p is the conformal ( "Liouville" ) factor of the worldsheet 
metric, is discretized as J VxVtp — > f Y[ v dx v . 

The planar diagram expansion is done with respect to the trivial (non-physical) vac- 
uum, M = 0. The propagator of the matrix field M with twisted boundary conditions 
is 

oo 

{Mf{x)M l j {x')) = G%(x,x') = e-l x - x ' +27r - Rm l(O m ) lJ (0- m ) fcZ . (4.4) 
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To each link I = (vv') of the planar graph is associated a propagator 



which is given by the m-th term in the expansion (|4.4|) . This gives rise to an integer- valued 
local field m = m< w /> on the graph. The weight of each planar graph (see Fig. 2) is given 
by 

tvfl w f ( \ 
W(E) = iV 2 - 2 V n^v" ex P - H \x v -x v/ + 2nRm <vv/> \\, (4.5) 

/££ V <vv>> J 

where n v is the total number of vertices in the graph; the integer field Wf is related to the 
field m< vv >> by 

Wf = rri£. (4-6) 
tedf 

where the sum is taken along the links on the boundary of the face /. 




Fig. 2: A discretized worldsheet. 



Geometrically, Wf is the winding number of the boundary df of the face / on the 
discretized worldsheet around the circle. By definition this is the vorticity associated with 
this face. Following standard arguments |36[] , the sum over the field me can be split into a 
sum over the vorticity field Wf and its gradient (vorticity-free) part.lll 

13 The latter sum, corresponding to m vv i = m v — m v / for some m v , m v i E Z, extends the 
integration over x v , which was originally restricted to the interval [0, 2ttR], to the whole real axis. 
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The dependence of a given planar graph on Q is through the factor Y[f g £ tvQ w f ; hence 
the twist matrix can be used to introduce fugacities of vortices with different vorticityB 
To study the effects of vortices on the worldsheet, one can either analyze the dependence 
of the partition sum ( }4.1| ) on Q (recall that in fact it depends only on the eigenvalues of 
Q), or integrate with respect to the twist matrix with an appropriate measure. We will 
use the second alternative, which is technically preferable. 

We start by showing that the singlet sector of the model is obtained by integrating 
over Q with the standard left/right invariant (Haar) measure [dQ] normalized so that 
f UN [dQ] = 1. An elegant way to evaluate the integral f UN [dQ] trO ni ■ ■ -trO nfc for ri\ + 
• • • + rik << N was proposed in f37j. It uses a Fock space representation of the moments 
of the unitary matrix in terms of bosonic oscillator modes. Replacing trO n — > a- n + a n , 
where {a n } and {d n } satisfy the canonical commutation relations [0^,0^] = [a m ,a n ] = 
m5 m+nj o, [ot-mi 0( n ] = 0, the above integral can be thought of as the expectation value with 
respect to the Fock vacuum a n \0) = a n \0) = 0. This leads to: 

! [dQ] tvQ n = N6 n , , 

(4.7) 

/ [dQ] trO m trO n = \m\5 m+nfi , etc. 

J U(N) 

Using these equations it becomes clear that the sum over all Feynman graphs (|4.5|) with 



a given set of vortices w f corresponds in the continuum language to the string worldsheet 
path integral with an insertion of the corresponding winding modes, ( [3.2| ): 



txQ n <-> T nR . (4.8) 



The 5-functions in eq. ( [4.7| ) impose winding number conservation. 
We would like to analyze the behavior of the partition sum 



[dQ] Z N (Q) (4.9) 

in the double scaling limit. By using ( ^.5| ) and the first line of ( ^.7| ) we see that each 
Feynman diagram with w/ = for all faces contributes a term that goes like N 2 ~ 2h (where 
h is the genus of the discretized Riemann surface). The continuum limit is achieved in 
the matrix model by tuning the coupling g in ( ^4.2| ) to a critical point g — > g c where large 
Feynman diagrams dominate the path integral, and considering the double scaling limit 
iV — > 00, g — » g c with N(g c — g) held fixed.lll The sum over diagrams corresponding to 
large worldsheets with Wf = scales like [(g c — g)N] 2 2h ~ /j, 2 ~ 2h in this limit. 



14 This generalization of the arguments of [27] was suggested by P. Zinn- Justin. 

15 up to certain logarithmic corrections, which are not essential here. 
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Feynman diagrams (|4.5|) with Wf 7^ do not contribute to the partition function in 
the double scaling (continuum) limit. In addition to the factor fx 2 ~ 2h discussed above, 
each insertion of a vortex is suppressed by 1/N, and the sum over graphs gives rise to a 
factor (g — g c )~ a where 1 — a = \uif\R/2 is the KPZ scaling dimension of the winding 
mode T Wf R (|3.2|). Since a < 1, the contributions of surfaces with vortices go to zero in the 
double scaling limit. 

Therefore, we conclude that the partition sum ( f4.9|) corresponds in the continuum 
limit to a sum over discretized worldsheets with no vortices, as claimed above. 

4-2. Introducing vortex couplings in the matrix model 

To enhance the contribution of vortices in the double scaling limit we replace the 
left-right invariant Haar measure [dfl] on the U(N) group manifold by 

[Dn) {x} = [DQ] exp j AntrO n ] (4.10) 

which is invariant only with respect to internal automorphisms O — > U~ 1 QU. The corre- 
sponding partition function is 

Z N [X] = j [DQ] {X} Z N {Q). (4.11) 



Applying the integration rules (|4.7| ) we see that in addition to the contractions discussed 



above, which vanish in the double scaling limit, we now have new terms where we use the 
second line of (|4.7| ) to contract a trO n from the exponent in ( |4. 10|) with one of the factors 



trO w/ in a Feynman graph ( f4.5| ). This leads to a factor of X Wf for each face with vorticity 
Wf. Therefore, the coupling constants A n can be thought of as fugacities of vortices with 
vorticity n (n = ±1, ±2, ...). 

One can also use the rules (|4.7| ) to contract the traces in the exponent in (|4.10| ) among 



themselves. This gives rise to a multiplicative factor exp J^ n>1 ?iA n A_ n in the partition 
function, which is large but "non-universal" in the language of non-critical string theory. 
To summarize, the free energy 

F N [\]=logZ N [\] (4.12) 

is given in the double scaling limit by a sum over all discretized worldsheets containing 
vortices, where the chemical potential of vortices with vorticity n is A n (up to a multiplica- 
tive factor which we will implicitly determine below). We conclude that the matrix model 
( |4.11| ) gives a discretization of the string theory ( |3.19| ). 
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4-3. Integrating over the twist angles in the double scaling limit 



The twisted partition function ( J4.ll ) can be calculated in the double scaling limit, 



where the matrix model simplifies drastically. As is familiar from the usual treatment of 
two dimensional string theory, in this limit the dynamics takes place near the top of the 
matrix potential V(M) (|4.2|) , and is dominated by large planar diagrams. Thus, to study 
the theory in the double scaling limit we can focus on the vicinity of the quadratic maximum 
of ([D| and replace the potential by V(M) = -M 2 /2 (after shifting M -> M - y/N/g). 
The behavior of the potential far from the top enters as a UV cutoff, which is usually 
taken to be a wall placed at a distance A ~ y/N / g from the top [|ioyTll , p!2| , p!3|Jl4l , p!5| ,|l( 
The twisted partition function (|4.1|) now reads 



Z N {Q)= [ VM e tr C Rd <-^ M ^ M2 ), (4.13) 

where the functional measure T>M is regularized as discussed above. 

For generic twist matrix Q we actually do not need to regularize the measure because 
the matrix integral is convergent. In the compactified theory, the cutoff has no effect when 
the compactification radius is smaller than half the self-dual radius (R < 1/2 in our units); 
we review this argument of [[27| in Appendix B. The Gaussian integration gives then a 
simple determinant depending only on the eigenvalues Zi,...,Zn of the twist matrix f2, 
which are phases, Zj = exp(z6>j): 



N 



where 

q = e 2l7lR . (4.15) 



Note that the partition sum ( |4.14[ ) depends on R only through q ( [4.15| ), and is also sym- 



metric under q — > the integral over the twist matrix (|4.11|) preserves this property. 
Naively, this means that R can be restricted to the interval < R < |. Later we will 
transform to the grand canonical ensemble, in which we will see that the symmetry is 
broken due to the appearance of certain divergences. 

The next step in evaluating the partition sum is the integration over Q in Q4.11D . Due 



to the symmetry of the measure with respect to internal automorphisms, the volume of 
the diagonal U(N) group factors out and the Q integral reduces to an integral over its 
eigenvalues zi, ...zn: 

[ [DQ] {X} Z N (Q) = ^f[f^ e 2u(zk) \A(z)\ 2 Z N (Q) (4.16) 
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where A(z) is the Vandermonde determinant (see e.g. Appendix A) and we introduced the 
potential 

u{z) = \Y,Kz n . (4.17) 

n 

Plugging ( [4. 14|) in (|4.16|) , we find that 

1 f N dz e 2u ( z k) N z . _ z ., 

ZN ^=MfU^t k (91/2-3-1/2) II 1/2 ' _ q -l/2 ■ ( 4 - 18 ) 
fe=l " jjtj' 

Another useful representation of Zn, obtained by using the Cauchy identity 

-dtf-U (4.19) 



Yli,j( a i- b i) (id)\ai-b 3/ 
is 

Zn[X] = TT / — det f eXP St )+M f ) • (4.20) 

A natural way to avoid the ambiguity in the contour integration in ( [4.20|) is to add a small 
imaginary part to the radius R and integrate over the {zk} along the unit circle.0 

The partition function can be calculated by residues for any N and but the formulae 
get more and more complicated when iV grows. Furthermore, the analytical structure of 
Zn makes it difficult to formulate a well defined 1/N expansion. 

4-4- The grand canonical partition function as a r-function of the Toda chain hierarchy 

The abovementioned problems can be overcome by considering instead of ( ^.18| ) the 
grand canonical partition function, where the cosmological constant \x becomes the chem- 
ical potential. Below we will show that the latter is a r-function of the infinite Toda chain 
hierarchy [[5S| and as such it satisfies the Toda equation (with appropriate boundary con- 
ditions), which plays the same role as the Painleve equation of the KdV hierarchy in the 
c = string theory. It will be convenient to rescale the vortex couplings and define new 
couplings 



A 



tn ~ q n/2_ q -n/2i ( 421 ) 



which are commonly used in the literature on r-functions. 
The grand canonical partition function is defined by 



Z^[\] = Y J ^ RllN Z n[\]. (4.22) 



7V=0 



16 Note that with this prescription we do not obtain a doubling of the degrees of freedom, 
typical for matrix models with symmetric potentials. 
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Using ( 4.20 ), one can show that the grand canonical partition function (|4.22|) can be 



represented as a Fredholm determinant 

Z^ = Det(l + e^ 2nR K), (4.23) 
where the integral operator K is defined as 

/J l u(z)+u(z') 
Ztci q 1 / z z — q L ' £ z 

In Appendix C, the partition function ( }4.23| ) is identified as a solitonic r-function of the 
Toda chain hierarchy.0 Explicitly the r-function with charge / is defined as 

oo 

n [t] = e -^ ntnt - e 27rRt *) N Z N [t] = e-^n nt ^Z„ +il [t] (4.24) 

N=0 

The grand canonical partition function ( [4.22]) is obtained by taking / = 0. The r-function 
with charge / can be obtained from ( ]4.22|) by replacing fj,—>fj, + il. 

4-5. Toda equation for the free energy 

The r functions r\ satisfy a hierarchy of Toda lattice equations which are differential 
equations in {t n } and finite difference equations in I. This infinite set of equations is 
generated by the corresponding Hirota identity, and reflects the Gl (oo) symmetry of the 
system. This symmetry can be exhibited by representing the system in terms of chiral free 
fermions. The lowest equation from the Toda chain hierarchy is 

nid+d-n) - (d+n^d-n) + n+tn-t = 0, (4.25) 

where we denoted d± = d/d\± and A± = t±i = ± 2i g^r^ ■ 

As mentioned above, the r-function depends on the cosmological constant \i and the 
charge / via the combination fx + il 

n(ji) = ro(jj, + il). (4.26) 

We will restrict ourselves to the case ( |3.20| ) (note that the constants A± are related to A±i 
by flPID ), for which 

r (^)=e- A + A -^(A). (4.27) 



17 Usually a solitonic r-function is denned by a finite sum instead of an integral. Our partition 
function is a particular case of the continuum limit (infinitely many solitons) of the soliton solution 



referred to as "general solution" in ref. [39j. 
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Winding number conservation implies that the r-function depends only on the product 
A 2 = A+A_, so that 

1 



d+d- = ^dxXdx, A = v^AT. 
The free energy fx) = logi? M (A) of the matrix model satisfies 

d+d.J^A, n) + exp \i + i) + F{\ H~i)- 2J r (A, //)] = 1. (4.28) 

In the double scaling limit .F(A, fx) is assumed to be a smooth function of \i. Thus we can 
rewrite (|4.28[) as 



d+d-TiX, \x) + exp I —4 sin 



\2dfi) 



F(\,fi)) =1. 



(4.29) 



Differentiating twice w.r.t. we get an equation for the susceptibility x(A, fj,) = 9 2 F(A, /x): 



d+d-x(\fi) + d 2 exp 



( sin(ia M; 



0. 



(4.30) 



The operator in the exponent is to be understood as a power series in d 



sin(ia M; 



i9 
2 W M 



-2n 



d 2 , n . 



n=l 



(2n+l)! M 



(4.31) 



Note that the differential equations (|4.29| ), ( |4.30|) do not depend on R explicitly. However, 
to solve them we need to supply boundary conditions. A convenient choice is to specify the 
free energy T at A = 0, and use the differential equations to solve for the A dependence. 
Since A = corresponds to the standard c = 1 matrix model, = T{\, /u)|a=o is 

known. We next review its form. 



4-6. Boundary conditions 

The c = 1 string theory without vortices corresponds to u(Q) 
integral in (|4.2Q|) can be evaluated by residues and gives 



Z 



N 



N 2 /2 



N 1 

n=l 



(see ( CT) ). The 



(4.32) 



This formula makes sense for the usual matrix harmonic oscillator, but for the inverted 
oscillator it requires an analytic continuation. Fortunately, the corresponding grand canon- 
ical partition function ( 4.22|) in the absence of vortices (all t m = 0) can be expressed as 
a partition function of fermions in terms of the density p(E) of the energy levels of the 
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inverted harmonic oscillator (see [ES|). The (universal part of the) grand canonical free 
energy of the c = 1 string theory compactified on a circle of radius R has a well defined 
(asymptotic) \j \i expansion: 

^) = logZ M [0] 

= - [°° dEp(E) log (l + e -**R(v-E)\ 

2 J -OO ^ ' 

= J_ f°° dEY* logfl + e~ 2 ^- E h (4 ' 33) 

fc=2 

where we introduced the polynomials in 

„. * f \\ 2n (2 2 ( fc -^) -2)(2 2n -2)|B 2(fc _ n) B 2 J 
MR) = (2fc - 3)! 2" M g (- J plLnlHW <4 ' 34) 

and _B m are Bernoulli numbers. The partition function (|4.33| ) has a T-duality symmetry, 
R — ► \i — > i?//. This symmetry is broken in the presence of vortices (i.e. for A 7^ 0). 

^.7. Canonical versus Grand Canonical partition sum 

The sum over worldsheets in the continuum string theory is described in the matrix 
model by the 1/N expansion of the canonical free energy jF^r = log^^v. The formalism 
described in the previous subsections allows one to compute the grand canonical free energy 
J-(lj) = logZ M , which is related to JF/v via 

exp T N = I — exp[2irRfiN + (4.35) 

J 27TZ 

where the integration contour encircles the point e~ 2nR ^ = 0. 

As is well known [^4| , |l5|Jl^| , the 1/N expansion of J-'nW can be rearranged as a 1/no 
expansion, with (i Q defined by (jlq log(/x /A) = N(g c — g). The grand canonical free energy 
JF(A = 0, /i) has al/jti expansion given by Q4.33Q . Similarly, in section 5 we will show that 
JF(A 7/ u) has a l/\i expansion which is obtained from ( |4.29| ). One can show that in the 
double scaling limit, the 1/no expansion of J-'n coincides with the expansion of J r (p), 
up to a flip of the sig of the leading (spherical) contribution to JF (see |f27| , revised 
version, section 6). 

In the next section we will construct the l/\x expansion of the grand canonical free 
energy, and use the above discussion to identify it with the genus expansion of the partition 
sum of the continuum theory. 
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This sign flip is a standard feature of the Legendre transform. 
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5. Some results on the genus expansion of the free energy 

The purpose of this section is to solve the differential equations (|4.29|) , ( |4.30|) with the 
boundary conditions ( |4. 33| ) . We have not found an exact solution of these equations, but 
as we will see below, one can solve them iteratively, in a genus expansion, 



r{\fj,) = Y t r h {\ t fj). (5.1) 

The genus h partition sum Th must, by KPZ-DDK scaling, have the formll! 



T h {\ fi) = S h , ^ 2 log^- 6 -^(R + i) log (jl + fi 2 - 2h A h (z), (5.2) 



where z is defined in (|3.9|). The boundary conditions (|4.33|) determine Ah[z) at z = 0. 



Plugging the ansatz ( |5T2| ) into Q4.29| ), ( |5.1|) one can solve for Ah(z). It is remarkable that 
the Toda equations allow such an ansatz at all; this is due to their conformal symmetry. 

At tree level one finds an ODE for Aq(z) which together with the boundary conditions 
fixes it uniquely. At genus h one finds a differential equation that mixes Ah(z) with Ah>(z) 
with h' < h. These equations can be used to determine all Ah iteratively given Ah(0). 

We will actually perform the calculations by transforming from the variables (/i, z) to 
(A, y), where y is defined in ( |3.17|) . In these variables the genus expansion is an expansion 



in inverse powers of X 2 - R , and the boundary conditions are specified at y = oo. The genus 
h free energy in ( |5.1|) takes the form 



(5.3) 



^i(A,/i) = ClogA + / 1 (y) 
T h {\ n) = (vX^) 2 - 2h fh(y) (h > 2), 
where, as in (|3.17| ), 



The genus expansion corresponds to fixed y with fi, A — > oo. The logarithmic term in F\ 
is a zero mode of the Toda equation ( f4.28| ). The constant C will be determined later. 

We note in passing that C is a special case of a large class of zero modes of the Toda 
equation ( [I.29|) , whose general form is 

oo 

A^=^(C n + J D n logA) e~ 2n ^. 

n=0 



19 The fact that the r.h.s. of ( [Op is one, leads to an additive contribution A+A- to T . This 
contribution is non-universal and will be dropped below. 



24 



For n^O these terms are non-perturbative in the genus expansion. The {C n } are nothing 
but the non-perturbative terms in the standard c = 1 string theory (|4.33|) . Their presence 
is related [§(J to the ~ (2/i)! divergence of the coefficients fh in the perturbative expansion 
( |5.3| ). The status of non-perturbative terms in two dimensional string theory is generally 
unclear. Here we are interested in the perturbative expansion of the free energy and will 
therefore drop all exponential corrections. 

Differentiating (|5.3| ) twice w.r.t fx we find the susceptibility, 



X 



(A, iA = ^2xh(x, n) 



Xo(A,/*) = log (Av^l) +X (y), ( 5 - 4 ) 

-2h 



X/ l (A,M)= (r/A^J X h (y) (h>l), 
with fh and Xh related as 

d 2 y f h {y)=X h {y). (5.5) 

5.1. Partition function on the sphere 

Plugging ( |5.4j ) into ( |4.30|) and retaining in the spherical approximation only the h = 
term (and only the leading term in ( |4 . 3 1|) ) we obtain the following ODE defining the string 
susceptibility X (y) 

a(yd y ) 2 Xo + d*e- Xo =0 (5.6) 
where a = (R — l)/(R — 2) 2 . The solution is 

y = e ~Ti x ° - e ^r x ^ (5.7) 



reproducing Moore's result ( |3. 1 1|) — ( |3.12| ). This expression can be checked by substitution 



or found by a direct solution of (|0|) (see Appendix D). The integration constant is fixed 
by comparing to the leading term in ( f4.33| ). 

Using ( |5.5|) and integrating Xq twice we obtain the following parametrization of the 
partition function on the sphere: 



Jo(A, n) = -\n 2 log(AV^T) - X ) 



3 R ofl-w 3„ _2 Xn R 2 -R + l _ x 



+ 2 A 4/(2-i?) -_^_ e -2^X + - Re ~iX 



AR-l 4 R-l 

with X (y) defined from ( |5.7| ). 



-e 



(5.8) 
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Consider the limits y — > 0, oo in (|5.8|) . In the limit y — > oo (fi — > oo or A — > 0) 
reproduces the known asymptotics of the c = 1 string theory perturbed by vortices: 

T Q ~ -^^ 2 log^-AV fl + 0(AV fl " 2 ), /u^oo (5.9) 

whereas in the Sine Liouville limit A — > oo with ,u fixed we obtain the asymptotics 

F ~ -AA 4 /( 2 "^ - S/xA 2 /^"^ + 0(// 2 ), A ^ oo. (5.10) 

where A = \{2 - R) 2 (R - 1) R /( 2 ~ R ) and B = (R - 1) V(2-fl)i|g£ § Note that the 
dependence of the partition sum on A is in agreement with the KPZ-DDK analysis in 
Sine-Liouville theory. In particular, for the black hole radius R = 3/2, we find JF ~ A 8 , 
in agreement with (|2.16|) . It is interesting that the partition sum on the sphere goes like 
an even integer power of A. This is a non-universal contribution to the partition sum of 
Sine Liouville theory, which corresponds to condensation of vortices. Recall that A is the 
fugacity of vortices, and a finite integer power corresponds to a contribution of a finite 
number of vortices and antivortices. Hence it should be dropped from the physical answer. 
As we shall see in section 6, this has an important physical consequence: the free energy 
of the two dimensional black hole vanishes in the tree approximation, in agreement with 



the results of [41 



Not also that the coefficient A is complex for R < 1, in agreement with our earlier 
analysis - the imaginary part for R < 1 signals the instability of Sine-Liouville (with real 
A) in this range. 

5.2. Linear equations for any genus 

Substituting ( |5.3| ) into ( [4.30| ) we obtain the following ODE for fh 

a(yd y + 2h- 2) 2 f h - e~ x ^d 2 y f h = W(/ , • • • , f h _{). (5.11) 

Here X Q (y) is defined by fl5.7p and the function TC(fo, • • • , fh-i) is defined as 

7~L(fo, - ", fh-i) = 

(v^TA)-^exp f-4sin 2 (|a M ) ^{^R^lX)-^^ f, 



k=0 



(5.12) 



o 



where by [• • -]o we denoted the constant term in the limit A, \i — > oo with y fixed. 

We see that the equations for {fh} form a triangular system of linear second order 
ODE's with y dependent coefficients. We do not know how to solve them in general; in the 
next subsection we will describe the solution for the case of genus one (the torus), where 
the equation simplifies. 



20 Note that a more natural parameter would be A = \/R— 1A; indeed, loking back at 

|3.9| ) we see that the dependence of the partition sum on A is through A. In terms of A one finds 
l 

4a' 



that^o^-^rA 4 ^ 2 -^. 



26 



5. 3. Partition function on the torus 

Eqs. (|5.11|) and (|5.12|) lead to the following first order linear inhomogeneous ODE for 
the torus partition function fi(y) 



W - e" Xo )/r + ayf{ = -i^e"^. 



(5.13) 



Note that ( |5.13|) is a first order differential equation for f[. Integrating this equation 
and using the boundary condition at large fi (|4.33| ) to fix the integration constant (see 
Appendix D) we obtain (up to a non-universal constant) 

R + R- 1 



24 



R Xq 



2-R 



1 

24 



log ( 1 - (R - l)e^ x ° 



(5.14) 

with X (y) defined by ( |5.7| ). In the limit /i — ► oo with A fixed we reproduce from ( |5.14| ) 
the known asymptotics of the c — 1 string theory perturbed by vortices: 



R~\~ R , R — 1,9 R 9 

log^+ — X jl 



+ ■ 



(fl -> oo) 



24 "~ or " ' 24R 

whereas in the Sine-Liouville limit A — ► oo and \x fixed we obtain the asymptotics 



(5.15) 



R + R' 1 
12(2 - R) 



log(XVR^l) 



R-l 



24(2 -R) ((R-l)A 2 ) 



l/(2-i?) 



+ 



(A -> oo). (5.16) 



The expansion ( |5.15| ) of the solution ( p.!4| ) can be compared with the explicit ex- 
pressions for genus one tachyon correlation functions in c = 1 string theory obtained in 



transformation R = 1/R, fx = Rfx. After this redefinition, the results of |24[ read: 

1 n\ 



24|j . Since [33J] studied momentum modes, to compare one has to perform a T-duality 

(5.17) 

r(r>o. - mi 

9n(R) 



where 



(T%7?R)h=i,R = ~ ^ [(i - R)n R ~ 2 ] n ((-) n fn(R) + g n R 2 ) 

T(n(2-R)) 



and 



T(n(l-R) + 1) 

f 1 (R)=R 2 -R-l, 

f 2 (R) = 3R 3 - 8R 2 + 3R + 3, 

f 3 {R) = 17R 4 - 72R 3 + 90R 2 - 17 R - 20, 

etc. 



For low n one finds 

(r + r_)^ R2 - 1 



24R 



R-2 



1 



12R 



1 - R) 2 (R S - 5R 2 +3R + 3) /i 



2R-4 



(5.18) 



This is exactly what one obtains by expanding the solution ( |5.14|) . 
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5.4- Comparison to continuum calculations 

Most of the matrix model results described above are difficult to derive using contin- 
uum methods. In this subsection we will discuss two cases in which the matrix model pre- 
dictions can actually be obtained directly in the continuum formulation; the non- vanishing 
of the partition sum on the sphere and the precise value of the partition sum on the torus. 

Consider first the genus zero case. In the matrix model we found a non-zero partition 
sum given by ( |5.10|) (with R = 3/2 for the Euclidean black hole case). In the continuum 
analysis on the cigarll this can be reproduced as follows. 

In string theory the partition sum on the sphere is usually said to vanish, due to 
the volume of the Conformal Killing Group (CKG) of the sphere, SL(2, C). If the target 
space is non-compact, the partition sum is actually proportional to V/vol(S r L(2, C)) where 
V is the divergent volume of spacetime. Thus, at first sight the partition sum is oo/oo, 
i.e. it seems to be ill defined. However, in most situations one is actually interested in 
the partition sum per unit volume. E.g. if the vacuum is translationally invariant in the 
non-compact directions, the partition sum per unit volume is the Lagrangian density in 
this vacuum (the classical cosmological constant) , and it vanishes due to the volume of the 
CKG. 

In the Euclidean cigar background /U(l), the above discussion has to be reexam- 
ined. There is no reason to divide by the volume of the cigar since the background is not 
translationally invariant in 0, and in any case for comparison with the matrix model we 
are interested in the total partition sum and not the partition sum per unit volume. The 
ratio of the volume of the cigar to the volume of the CKG is finite in this case. The volume 
of H$ contains precisely the same kind of divergence as that of the CKG. Since the volume 
of the U(l) in (|2.1|) is finite, we conclude that the partition sum of string theory in the 
cigar background is non-zero (and proportional to the mass of the black hole M). We will 
not attempt to compute its precise value. 

Moving on to the torus, we would like to derive the result (|5.16|) directly in the 



continuum formalism. We will describe the calculation in the Sine-Liouville language, to 
study the R dependence. Similar considerations can be used to perform the calculation on 
the cigar. 

The torus partition sum in the background ( |2.7| ) is an example of a bulk amplitude 
in the language of section 2. By performing the integral over the zero mode of (p as in eq. 
( |2.21| ) and using ( |2.23| ) one finds 



1 , A f d t 



'torus — 2 — R A ' ^ 2 



log T / -^-Zi(t,t), (5.19) 
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A similar analysis can be performed directly in Sine-Liouville theory. 



28 



where Zi(r, f) is the torus partition sum of the theory with A = 0, i.e. string theory on an 
infinite cylinder of radius R. Z\ can be computed by using free field techniques; the linear 
dilaton term in (|2.7|) can be ignored in this case since the curvature of the torus vanishes. 
Equation ( |5.19| ) has a simple interpretation: the bulk of the contribution to the torus 



partition sum comes from the region far from the Sine-Liouville wall. It can be computed 
by cutting off the infinite region very far from the wall, by restricting to (j> < (puv- As 
indicated by the notation, this can be thought of as adding a UV cutoff to the theory. 
The Sine-Liouville wall acts as an effective IR cutoff, restricting 

Xe -(2-R)<p < ^ or ^ > _}_\ og \ (5.20) 

2 — R 



Thus, the length of the cylinder is 



'UV 



2-R 



log A. (5.21) 



The prefactor of the modular intgeral in (|5.19|) is hence nothing but the (universal part 



of the) length of the cylinder, and the integral computes the free energy of perturbative 
string modes (tachyons) living on the cylinder. 

This makes it clear that the coefficient of — log A in ( |5.19| ) must be the same as 
the coefficient of — ^ log [i in the torus partition sum in standard c = 1 string theory, which 
has been computed by performing the integral ( |5.19|) (see e.g. fl5|l ). As we see in ((L33|), 



this coefficient is [R + 4)/12, in agreement with ( |5.16| ). 



6. Thermodynamics of two dimensional string theory 

Our results from the previous sections can be used to study the thermodynamics of two 
dimensional string theory. We start by reviewing the low temperature thermodynamics 
(corresponding to R > 2 in the notations of the previous sections), and then describe 
the physics for temperatures slightly below and slightly above the Hagedorn temperature, 
which corresponds to R = 3/2. We also comment on the behavior in the intermediate 
region 3/2 < R < 2. 

6.1. Low temperature thermodynamics, R > 2 

The canonical partition sum Z(f3) is obtained by computing the string path integral 
with Euclidean time compactified on a circle of circumference f3 = 2-kR. The precise 
relation is 

log Z((3) = -PF{0) = Z string (R), (6.1) 
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where F is the free energy and Z str [ ng is the string partition sum, which is given by a sum 
over connected Riemann surfaces. For two dimensional string theory, below the Kosterlitz- 
Thouless temperature (i.e. for R > 2 in our units) one can ignore winding modes, i.e. set 
t n = in (|3.19| ), and the free energy (|6.1|) is given by the familiar formula (|4.33|) 

Z stliag (R) = -| n 2 log | - ±(R + i) log £ + • • • . (6.2) 

We have written the first two terms in the 1/// (or genus) expansion. A is a large scale 
A >> fx corresponding to a UV cutoff. 

Comparing to (|6.1| ), we see that the free energy F(i3) is given by 



The temperature independent terms in (|6.3| ) can be attributed to vacuum energy. Indeed, 
the energy 

goes as R — > oo to the ground state energy 

and it is natural to subtract it from both E and F. This leads to 



1 1 , V 

E = -^24Ri l ° g A + 

F = H log — + 



(6.6) 



We can now also use the thermodynamic relation 

—(3F = S-f3E (6.7) 

to determine the entropy, 

S = /3(E-F) = -^log^- + .... (6.8) 

Eqs. (|6.6|) , (|6.8|) have an obvious interpretation: they correspond to thermodynamics of a 
single massless 1 + 1 dimensional field which lives in a spatial volume 

V L = -log^. (6.9) 
Indeed, solving fl6.6| ) for R in terms of E, Vl, and substituting in (|6.8| ), one finds 



S-^ (0,0, 

which is the standard result for a single massless scalar field in 1 + 1 dimensions. The scalar 
field in question is of course the massless "tachyon" of 1 + 1 dimensional string theory, 
the only perturbative field theoretic degree of freedom. Thus, we conclude that the low 
temperature thermodynamics corresponds to a gas of tachyons. 
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6.2. Near- Hag edorn thermodynamics, A: R> 3/2 

As discussed in the introduction, the high energy density of states of two dimen- 
sional string theory has the Hagedorn form ( |1.2[ ). Thus, it is interesting to consider the 
thermodynamics in the vicinity of the Hagedorn temperature Tjj = 1//3h- As we show 
below, f3 h = 2ttRh, with Rjj = 3/2. The thermodynamics in the vicinity of the Hagedorn 
temperature is quite sensitive to the value of the constant si in (|1.2f) . We will see below 
that 

-2 < si < -1. (6.11) 

For now we will take (|6.1l| ) for granted and study its consequences for the thermodynamics 
at temperatures slightly below Tjj. In the next subsection we discuss the situation for 
temperatures slightly above Tjj, and in the process determine both Tjj and si. 
Consider the canonical partition sum 

Z({3) = / dMp(M) e~ m . (6.12) 



The high energy behavior of p(M), given by (.L2), implies that the integral converges for 



T < Th, but the exponential Boltzmann suppression disappears as (3 — > (3h- Naively, 



one might expect that in this limit the integral (|6.12| ) becomes dominated by high energy 
states, so that 

Z((3) ~ ((3 - (3 H )- s i-\ (6.13) 

but this is not quite the case. The contribution of the high energy part of the spectrum to 
Z is given in (|6.13| ). The rest of the states give an additive contribution Z\ ow ((3) which is 
analytic as (3 — > (3h and approaches a non-zero constant in this limit. The partition sum 
( |6.12| ) can thus be written for (3 slightly larger than (3h as 

Z{(3) = Zum(fi) + c{(3 - (3 H )~ S1 ~\ (6-14) 

where c is a constant. Since si + 1 is negative (see ( |6.11[ )), the contribution of the high 
energy states (the last term in (|6.14|) ) goes to zero, and the partition sum (|6.12| ) is in fact 
dominated by states with moderate energies. 

Naively, the Hagedorn temperature appears to be a limiting temperature, since the 
energy as measured in the canonical ensemble diverges: 

E = -^f~ {( 3-f3 H) -^. ( 6 .15) 

Note that —si — 2 is negative (see ( |6.11| )); thus the derivative of the second term on the 
r.h.s. of (|6.14j ) is much larger than that of the first. Hence, it appears that one needs to 
supply an infinite amount of energy to heat up the system to the Hagedorn temperature. 
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However, one has to be careful with this conclusion, since the whole notion of the 
canonical ensemble breaks down as one approaches the Hagedorn temperature from below, 
since the fluctuations around the mean energy are large: 

~(/?-ftr)' 1+1 - (6-16) 

In such situations the system should be studied in the micro canonical ensemble, where the 
energy is fixed, and the temperature is given by 



Using the entropy (|1.2| ) we find 



P = P* + £ + o(±). (6.18) 



M \M 2 

Since s\ is negative (see (|6.11| )), we find that at high energies the temperature is above 
the Hagedorn temperature. Clearly, this does not mean that temperatures below the 
Hagedorn one are not achievable in two dimensional string theory - as we reviewed in the 
previous subsection, the low temperature thermodynamics is perfectly sensible. Instead, 
( |6.18|) implies that the microcanonical thermodynamics slightly below Tjj corresponds to 
a finite energy M, obtained by solving ( |6.17|) . 

Thus, if we imagine a process where two dimensional string theory is gradually heated 
up by supplying energy to the system in the microcanonical ensemble, the temperature Th 
is reached at a finite value of the energy M. This suggests that the Hagedorn temperature 
is associated with a phase transition (as opposed to being a limiting temperature) . 

6.3. Near- Hagedorn thermodynamics, B: R < 3/2 

In the previous subsection we saw that just below the Hagedorn temperature the 
thermodynamics is dominated by states with moderate energies. Since two dimensional 
black holes should only dominate the thermodynamics at large energies, it is clear that the 
low temperature phase is not well described by perturbative string theory in the Euclidean 
black hole background ( |1.1| ), even as T — > Th- Instead, as we show in this subsection, string 
theory in the background ( |1 . 1| ) describes the thermodynamics just above the Hagedorn 
temperature. 

For very large M, the thermodynamics is obtained by studying classical string theory 
on the cigar. The mass M is related to the value of the string coupling at the tip of the 
cigar, while the temperature is determined by the asymptotic radius of the cigar. Since 
the latter is independent of the former, we conclude [IT] that 
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is constant and the entropy-energy relation is 



S = (3 H M (6.20) 

i.e. the density of states exhibits Hagedorn growth, p(M) ~ exp(/3#M) at large M. The 
free energy ( |6.7| ) should thus vanish in the classical approximation. 

Quantum string theory on the cigar leads to 1/M corrections to ( |6.19| ), ( |6.20|) , whose 
general form is (see (|1.2|) ) 



S=0 H M + 8l ]ogM + o(Jj 



M \M 2 



(6.21) 



where si, a are constants to be determined. 

For [3 < Ph, the integral (|6.12j ) defining the canonical partition sum diverges. At first 
sight this looks like a problem, but in fact it is standard in black hole theory. For example, 
black holes in flat spacetime generally have an entropy that grows like M a with a > 1, so 
the analog of ( |6.12| ) diverges for all temperatures. The Gibbons-Hawking prescription [42 



whose generalization to string theory we are studying here, defines the canonical partition 
sum by a formal saddle point evaluation of the integral ( |6.12| ) . The imaginary part of the 
resulting partition sum is associated with the thermodynamic instability of these black 
holes (their negative specific heat), and is related to the appearance of an unstable mode 
(a tachyon) in the Euclidean black hole background. 

We will see that something very similar happens here. Two dimensional string theory 
is thermodynamically unstable for T > Th- Nevertheless, just like in ||42|| , we can compute 
the free energy by a formal saddle point evaluation of ( |6.12| ). The Euclidean black hole 
background has a tachyon, whose condensation needs to be understood in order to study 
the high temperature phase of two dimensional string theory. 

It is not difficult to show that a saddle point evaluation of the integral (|6.12[) with the 



density of states ( |L2"| ) gives the result ( |6.13| ). The saddle point is located at M given by 
( |6.18| ), and assuming that si is indeed negative ( |6.11| ), as we will find soon, we see that for 
(3 slightly below (3jj it is located at a large positive value of M, which is consistent with 
our use of the asymptotic density of states (|1.2|) to find it. 

Since the partition sum is given by (|6.13|), the free energy is 



-f3F = -{ Sl + 1) logOS - M + 0(1). (6.22) 

The energy-temperature relation (in the canonical ensemble) is obtained from ( |6.22|) by 
using 

M = «m = |i±l. (6 . 23) 

d/3 f3 - 13 H 
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Note that it is different from the relation ( |6.18|) , which is valid in the micro canonical 
ensemble. The two are not the same here because of the large energy fluctuations in the 
canonical ensemble. 

Substituting (|6.23| ) in ( |6.22|) we conclude that 



-QF= (si + l)logM + 0(l). (6.24) 



Thus, by computing the string partition sum (|6.1| ) we can deduce the entropy-energy 
relation characterizing the black hole, ( |6.21| ). 

An immediate puzzle is that while we expect the tree level contribution to the free 
energy (the term that goes like M in ( |6.24[ )) to vanish, our explicit calculation in section 
5 gave in fact a finite answer for this coefficient. This apparent paradox is resolved by 
noting [4~I] that the thermodynamic relation ( |6.7| ) is modified in this case by the presence 



of a massless scalar field, the dilaton, that mediates a long range force. The black hole 
spacetime is in fact charged under the dilaton current J a = e^Vf,e^, and therefore the 
thermodynamic relation ( |6.7|) contains a correction that is proportional to the dilaton 
charge of spacetime. It was shown in [El] that this gives rise to a term in the free energy 



that goes like M (times a constant that is ambiguous in the gravity approximation), and 
should be interpreted as a chemical potential for the dilaton charge. 

Therefore, we believe that the correct procedure for studying the thermodynamics of 
two dimensional string theory is to subtract this term from F. Note that in the Sine- 
Liouvile description, this is a term that goes like A 8 (see (|5.10|) ). In general in non- 
critical string theory, such terms are due to physics far from the Liouville (or in this 
case Sine-Liouville) wall and are "non-universal" (as noted at the end of section 5.1). 
In our particular case one could alter the coefficient of this term by adding a constant 
(independent of to the eight point function (T^T^ R ) given in ( j3Jj|) ; such a constant is 
non-universal in Liouville theory (or the double scaled matrix model). 

As further justification for the validity of this procedure we discuss in the next section 
a generalization of two dimensional string theory to the fermionic case (with worldsheet 
supersymmetry but no spacetime supersymmetry) . The physics of this theory is very 
similar to that of the bosonic one, but the tree level partition sum vanishes. 

The one loop correction to the free energy discussed in section 5.3 (see (|5.16|) ) implies 

that 

si + 1 = -( J R + .R" 1 )/48 = -13/288. (6.25) 

which is in agreement with ( |6.11|) . 

The high temperature thermodynamics described by the black hole is unstable. As 
is clear from ( |6.23| ), ( |6.18| ), the specific heat is negative: as the energy increases, the 



temperature decreases towards Tjj- This is of course closely related to the fact that in 
the formal saddle point evaluation of the partition sum (|6.12| ), we are in fact expanding 
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around a maximum of f(M) = f3M—S(M) (rather than a minimum) . This thermodynamic 
instability is associated with the presence of a negative mode in the Euclidean path integral, 
i.e. one expects to find a tachyon in the cigar background ( |lTT|) . 

Because the instability of the thermodynamics is a one loop effect in the cigar back- 
ground, one actually expects to find a zero mode in the classical theory, which becomes a 
negative mode after one loop effects are included. There indeed is a very natural candidate 
for precisely such a mode which we will describe next. This issue was recently discussed 
in the context of LST in || , so we will be brief. 

It is known (see e.g. ||) that in the cigar geometry there are normalizable states 
localized near the tip, corresponding to the principal discrete series representations of 
SL(2), i.e. for 

\m\=j + l; 1 = 1,2,3, (6.26) 

One way to find these states is to study the correlation functions of the non-normalizable 
observables Vj; m ,m (see section 2); the states ( |6.26| ) give rise to poles in these correlation 



functions. The classical zero mode (or "massless state") referred to above corresponds t 

m = to = |; j = | - 2. (6.27) 

It is easy to check that this combination of j and m satisfies Aj^^ = Aj^^ = 1 (see 
(|2.3|) ), i.e. this is an on shell masless state (a zero mode). As discussed in 0, it is likely 
that including one loop effects turns it into a negative mode. It would be interesting to 
prove this by an explicit one loop calculation. This should be possible using the results of 
section 5. 

Thus, we conclude that our perturbative analysis of the black hole background de- 
scribes the high energy phase of two dimensional string theory. The thermodynamics is 
unstable, and there is a negative mode (or more precisely a classical zero mode that is 
conjectured to become unstable when one loop effects are taken into account). The high 
temperature phase presumably involves condensation of this mode and would be interesting 
to study using the techniques of sections 4, 5. 

6.4- The intermediate region 3/2 < R < 2 

So far we have discussed the thermodynamics below the Kosterlitz-Thouless temper- 
ature, i.e. for R > Rkt = 2, where it is dominated by a gas of perturbative string states 
(see ( |6T6| ) - ( |6.10| )), and at temperatures near the Hagedorn temperature. This leaves an 
unexplored intermediate region, 3/2 < R < 2. We will not study this region in detail here, 
but would like to make a few comments on it. 



22 Since the result is general, we are writing it as a function of k. For our case, as usual, k = 9/4. 
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As the temperature is raised above the Kosterlitz-Thouless one, the system is ex- 
pected to undergo a phase transition to a phase where the non-singlet degrees of freedom 



contribute to the thermodynamics ||25|1 . In this phase, the free energy and entropy are 
expected to receive tree level contributions (in contrast to (|6.6|), (|6.8| )), associated with the 
large density of non-singlet states. 

In the continuum description, one expects the winding modes ( |3.19|) to condense. For 
1 < R < 2, only the winding number one mode is relevant, and thus we expect t±± = A/2 
in (|3.19| ) to be the only non-zero couplings, in addition to \i. An important question for 
analyzing the thermodynamics is how does A depend on R. It is not clear to us how one can 
determineil A(i?), but it appears likely that it is such that the parameter z defined in ( |3~fl| ) 
varies smoothly between zero and some finite value, as R decreases from R = 2toR = 3/2 
(the Hagedorn temperature). This variation is related to the dependence of the energy 
on the temperature. As R — > 3/2 the energy should approach a finite value (as discussed 
above), while as -R — > 2 we should connect to the low temperature thermodynamics, for 
which z = 0. 

The qualitative picture proposed in the previous paragraph is in agreement with 
the phase diagram of the Sine-Gordon model (or, in statistical mechanics language, the 
Kosterlitz-Thouless phase diagram). In the region R < 2 the phase diagram is such that 
as the distance scale increases (i.e. as one goes from the UV to the IR), A grows and R 
decreases. String thermodynamics for 3/2 < R < 2 can be thought of as the Sine-Gordon 
model coupled to worldsheet gravity. It is well known that the qualitative structure of 
the RG in two dimensional QFT coupled to gravity is the same as in the theory without 
gravity, with \i playing the role of the RG (energy) scale, /i — > oo corresponds to the UV, 
while small \i is the IR region. Thus, in the Sine-Gordon model coupled to gravity one 
expects to find A = and R — R(/j), i.e. A = X(R). It is not difficult to see that the 
qualitative properties obtained from this picture are in agreement with the proposal in the 
previous paragraph. A more precise analysis is left for future work. 



7. Fermionic two dimensional string theory 

In this section we briefly describe a natural generalization of two dimensional string 
theory to the fermionic (type 0) case. Our main motivations for discussing this theory are 
the following. 



23 Naively, A, R are two independent parameters that specify superselection sectors in the 
theory. However, as we saw above and discuss further in section 8, it is important here to keep 
the UV cutoff A finite. This makes A, R fluctuating fields, and one must extremize with respect 



to them. See []43|] for a related discussion. 
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In section 6 we argued that in the bosonic 2d string, the tree level free energy that 
goes like the mass of the black hole, corresponds to a vacuum energy type contribution, 
and should be subtracted. We will see that most of the physical properties of the type 
theory are very similar to the bosonic one, but the tree level free energy vanishes. This 
provides further evidence for the validity of the procedure of section 6. 

The second reason for studying this generalization is that it provides a bridge to higher 
dimensional physics. In particular, the two dimensional worldsheet supersymmetric black 
hole plays an important role in the study of the thermodynamics of Little String Theory 
(see for a recent discussion). 

Two dimensional type string theory contains two worldsheet superfields, X, with 
central charge c = | and <fi, with c = | + 6Q 2 . The fact that the total central charge is 
c = 15 determines Q = y/2. The zero temperature physics of the theory is described in 
The perturbative spectrum includes in this case two massless scalars, a NS-NS field 
T(X), and a R-R field V(X). As in the bosonic case, there are no perturbative massive 
string states (for generic momenta). The bulk amplitudes of T and V are known; they are 
completely determined by an infinite dimensional symmetry algebra, which is very similar 
to that found in the bosonic string, and are themselves very similar to the bosonic ones. 

Our main interest here is in the high temperature thermodynamics of the model. Just 
like in the bosonic case, it should be dominated by states associated with 1 + 1 dimensional 
black holes. The black hole background coresponds in this case to an iV = 1 super conformal 
version of the coset (|2.1|). Algebraically, it can be thought of as follows. 



One starts with the N = 1 superconformal theory on , which contains an 
SL(2,R)l, x SL(2, R)r current algebra with total level k. The level k receives a con- 
tribution of k + 2 from a bosonic WZW on H$ , and —2 from three free (left and right 
moving) fermions in the adjoint of SL(2). The total central charge of the theory is 

3(^ + 2) 3 

c H = + -. (7.1) 

One then gauges a U(l) super-Kac-Moody algebra, thereby decreasing the central charge 
to 

_ 3(fc + 2) 

t'coset — k ' \ ) 

The requirement that the total central charge is equal to 15 determines 

11 

Q 2 " 2 

Much of the discussion of section 6 goes through here. The classical analysis gives^ 



k = ^ = -. (7.3) 



S = (3 H M; p H = 2ic~. (7.4) 



24 We work in the same units as in the bosonic theory. The Kosterlitz-Thouless radius is R — y/2 
in the fermionic case. 
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Therefore, the classical free energy ( |6.7|) should vanish. Unlike the bosonic case, here this 
can indeed be shown to be the case, as we briefly review next (see || for further discussion). 

As discussed in section 5, in the bosonic Euclidean cigar background, the spherical 
partition sum (which is related to the classical free energy via ( |6.1|) ) is in general non-zero. 
The volume of the cigar cancels the volume of the Conformal Killing Group of the sphere 
and leaves behind a finite residue. 

In the superstring one has to take into account fermionic zero modes. The CKG is 
replaced in this case by a supergroup, SCKG, but the associated fermionic zero modes are 
cancelled by fermionic zero modes of the N = 1 SCFT on the cigar. Thus, it appears that 
in the fermionic string, the partition sum in the background (|2.1[) is non-zero as well. 

This conclusion is incorrect because of the fact, mentioned at the end of section 2, that 
the SCFT on the cigar has in fact a larger symmetry, an accidental N = 2 superconformal 
symmetry. Thus, this SCFT has twice as many fermionic zero modes as one would guess, 
and integrating over them leads to a vanishing spherical partition sum . 

The vanishing of the free energy in this case is due to N = 2 worldsheet supersym- 
metry. It does not require spacetime supersymmetry, which in fact is absent in the type 
theory under consideration. It is also interesting to note that the free energy only vanishes 
at the Hagedorn temperature. This can be seen by using the duality of relating the 
Euclidean cigar SCFT to iV = 2 Liouville theory. As discussed in [Q for k < 1, which is the 
case here ( |7.3| ), the N = 2 Liouville is more weakly coupled and thus more appropriate. 

As one changes the radius of Euclidean time, the N = 2 Liouville superpotential 
changes to one that only preserves an iV = 1 subalgebra of the N = 2 superconformal 
symmetry. Thus, the arguments above no longer imply that the free energy vanishes, and 
one expects that it does not. This is in agreement with the thermodynamic picture that 
one expects: the free energy goes to zero as the energy goes to infinity, or (3 — > (3h- At 
finite energy the free energy is non-zero. 

One can also repeat for this case the calculation of the one loop partition sum in 
section 5.4, and the discussion of section 6, with very similar conclusions. To go beyond 
one loop probably requires finding a matrix model that describes the dynamics. This is 
an interesting open problem. 



8. Some open problems and future directions 

In this paper we discussed the thermodynamics of two dimensional string theory above 
the temperature corresponding to the Kosterlitz-Thouless phase transition. In the Eu- 
clidean time formulation of the problem, this regime is characterized by the condensation 
of winding modes around Euclidean time. We showed that the theory undergoes a finite 
temperature phase transition around the Hagedorn temperature Tjj. The high tempera- 
ture phase near T ~ Tjj is well described by the Euclidean two dimensional black hole. 
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Furthermore, by using the conjectured correspondence [21] of string theory in the black 
hole and Sine-Liouville backgrounds, we constructed a matrix model which provides a pow- 
erful tool for studying the two dimensional black hole to all orders in string perturbation 
theory. 

There are many open problems associated with our work, some of which are listed 
below. 

While the differential equation (|4.29|) in principle determines the partition sum ^(A, //) 
to all orders in string perturbation theory (as described in section 5), we have so far not 
been able to solve the resulting equations ( |5.11| ) beyond genus one. Furthermore, ( |4.29j ) 
is just the lowest equation in an infinite hierarchy which should allow the calculation of 
correlators of winding modes with arbitrary winding number at any order in g s . Some low 
genus correlation functions, like ( |2.5| ), have been computed in the continuum formalism, 
and it would be nice to reproduce them using the Hirota equations (see appendix C). These 
calculations remain an interesting open problem. 

There are observables in two dimensional string theory which seem to be outside of the 
framework of the Toda equations. Examples include modes carrying momentum around 
Euclidean time, and discrete states. It would be interesting to generalize the formalism to 
include these couplings. 

In the discussion of black hole thermodynamics in section 6, we assumed that the 
mass of the black hole is related to the string coupling at the tip of the cigar via the 
classical relation M ~ exp(— 2$ )- This relation is expected to be corrected quantum 
mechanically, and such corrections will influence the thermodynamic relations at two loop 
level and beyond. Thus, this issue will need to be addressed when the higher order results 
mentioned above become available. 

It would be interesting to understand better the thermodynamics of two dimensional 
string theory in the intermediate region 3/2<i?<2, as well as in the high temperature 
phase R < 3/2, and in particular to explore further the relation to the Kosterlitz-Thouless 
phase diagram mentioned in section 6. 

The FZZ correspondence ||21|| described in section 2 played an important role in our 
discussion, since it allowed us to replace the curved geometry of the cigar, which is difficult 
to realize in MQM, by the winding mode condensate which is easier to study in the matrix 
model. This correspondence remains very mysterious. It is not clear how the metric of the 
cigar is encoded in the Sine-Liouville potential. The phenomenon is probably much more 
general {e.g., as we mentioned, it can be extended to the fermionic string) and should be 
very interesting to understand better. 

String thermodynamics seems to suggest that two dimensional string theoryU has a 



Hagedorn spectrum of states at high energies, given by eq. (|1.2|) . This is surprising, since 



25 Note that here we mean the standard c = 1 string theory at zero temperature and with the 
scale set by the cosmological constant /x. 
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the perturbative spectrum consists of one field theoretic degree of freedom. Thus, it must 
be that most of the states described by (|1.2| ) are non-perturbative, and it would be very 
interesting to understand what they are .1^1 A natural proposal is that these states belong 
to the non-singlet sector of the matrix model. The degeneracy of non-singlet states is very 
large, but their masses are proportional to log(A//x) (with A a UV cutoff) |35| and thus 
diverge in the limit A — > oo. It appears that one can define two dimensional string theory 
in two different ways. 

If one sends A to infinity, all the non-singlet states are removed from the spectrum. The 
matrix dynamics reduces to that of just the singlets (i.e. the eigenvalues of M ( |4.1|) ), which 
in spacetime means that the massless "tachyon" is the only physical degree of freedom. 
The low temperature thermodynamics reviewed in the beginning of section 6 applies then 
for all temperatures, and all the vortex couplings t n in ( |3.19| ) are set to zero. We believe 
that the resulting theory does not admit two dimensional black holes (|1 . 1| ) and does not 
have a Hagedorn density of states. 

A second possibility is to keep the UV cutoff A large but fixed. Then the energy gap 
to the non-singlet states remains finite, and the physics is much richer. In particular, the 
model does contain black holes, and the corresponding Hagedorn spectrum of states (|1.2| ) 
should be reproduced by counting non-singlet states in MQM. 

It may seem puzzling that there is more than one way to define two dimensional string 
theory, but in fact this situation is standard in holographic theories that were studied in 
recent years. Consider, for example,ill string theory on AdS3. As is well known, three 
dimensional gravity with negative cosmological constant can be quantized as a Chern- 



Simons theory. According to the AdS-CFT correspondence [4"3|, this theory is equivalent 
to a two dimensional CFT with a large central charge (in the semiclassical limit), which 
in this particular case contains only the conformal block of the identity. The high energy 
density of states of this theory is much smaller than that given by (BTZ) black hole 
thermodynamics; thus this theory does not contain 2 + 1 dimensional black holes. 

An alternative definition of gravity on AdSs is obtained by studying string theory 
on AdSs and including all the perturbative and non-perturbative string excitations. This 
gives rise to a much richer theory, and in particular one can show that it does have black 



26 We note that a very similar situation occurs in weakly coupled Little String Theory in the 
double scaling limit defined in [||. The perturbative spectrum has a Hagedorn density of states 



of the form (1.2), but thermodynamics suggests that the full theory has a richer spectrum. In 
particular, (3h is predicted to be larger than the perturbative one. Thus, most of the states are 
non-perturbative. 

27 Similar comments can be made for d > 3, but the case d = 3 is special in that three 
dimensional gravity without matter exists as a quantum theory, and that the theory has an 
infinite dimensional symmetry, like two dimensional string theory. 
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hole states and their entropy can be computed microscopically, by studying the CFT dual 
to string theory on AdS^. This CFT can be formulated (roughly) as the extreme infrared 
limit of a matrix field theory; the bulk of the black hole entropy comes from off-diagonal 
components of the matrices. Moreover, just like in our case, the mass gap to the lowest 
lying black holes is non-perturbative in the string coupling. 

In two dimensional string theory, the above analogy suggests that we should keep A 
large but fixed in the double scaling limit. This is in fact very natural since A ~ N, so 
we are proposing that the non-perturbative theory has large but finite N. This is again 
similar to string theory on AdS space, where the analogs of N, [i are kept large but fixed.il 

The large entropy ( |1.2|) should then be due to non-singlet states. Not much is known 
about the spectra of such states at present, especially in the large N limit of interest here. 
One approach to this problem is to rewrite the partition function ( f4.11| ) as a sum of Gibbs 
partition functions over representations r of SU(N): 

Z N ((3,X) = Y,9r[X]Tr r e-^, (8.1) 

r 

where 

g r [X] = f [Dn} Xr (rt)exp ( ^ A n trfi n J , (8.2) 
(x r (fit) is the Weyl character in the representation r) and 

^ N 1 f) 2 1 1 T T -T r - 

Hr = Pr y r _ _ + L y (8 . 3 ) 

r f^ 1 2dx 2 k 2 ki 2 ^ (xi-Xj) 2 V ; 

k=l ^ J 

is the Hamiltonian in representation r in terms of the eigenvalues Xi of the matrix M with 
the inverted oscillator potential. This Hamiltonian is a matrix in the representation space: 
P r is a projector to this space and t[a are SU(N) generators on r (see |2^,^5j for details). 



These Hamiltonians define, in principle, in a natural way the time evolution of the system. 
Therefore to count the states in a given energy interval one can try to diagonalize them for 
representations corresponding to large Young tableaux. This appears to be difficult due to 



the Calogero type interaction of the eigenvalues with an SU(N) spin structure (see [ p6| , p7 ] 
for a discussion of the adjoint representation). 

If the non-singlet states indeed account for the entropy of black hole ( |1.2[ ) , it is natural 
to ask whether one can think of non-singlets as excitations of black holes directly. A 
possible direction for making the relation more precise is to note that the non-singlet states 
can roughly be thought of as open strings. In black hole physics it has been proposed in 



28 These analogs are the radius of curvature of the AdS space in string units, and in Planck 
units. 
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the past that the black hole entropy is due to open strings whose ends lie on the horizon 
of the black hole. It would be interesting to make this relation more precise. 

An interesting set of issues is associated with the Minkowski continuation of our 
results. The description of the black hole in terms of the Sine-Liouville model involves 
winding modes around Euclidean time and thus seems inherently Euclidean. Yet, many of 
the interesting questions about black hole physics can only be discussed in the Minkowski 
version. Examples are the unitarity of scattering off black holes and the understanding 
of the dynamics of black hole formation. In principle, the Hamilonians ( |8.3| ) provide a 
direct way for the analytical continuation to Minkowski space. Thus, to understand the 
formation and Minkowski evolution of black holes one needs a better understanding of the 
dynamics of the non-singlet states. This will have to be left for future work. However, we 
would like to comment on one point related to black hole formation in two dimensional 
string theory. 

It is sometimes stated in the old matrix model literature that the black holes ( p. . 1|) are 
not good toy models of higher dimensional black holes, since they cannot be formed in two 
dimensional string theory because of the large symmetry present in this theoryJH This 

symmetry (see [ P~5|JT^ ] ) is closely related to the interpretation of tachyon dynamics in 
terms of free fermions, with the tachyon describing ripples on the Fermi surface. The fact 
that the underlying fermions are free does not allow black holes to form in the evolution 
of tachyon pulses, even if they are large. 

Our discussion in this section clarifies the extent to which this statement is expected 
to be correct. In the theory with A — > oo, which corresponds to free fermions and describes 
the dynamics of tachyons as perturbations of the Fermi surface, one indeed does not expect 
to be able to make black holes since the theory does not contain black holes. However, 
in the theory with finite A, the free fermions are just one sector of the Hilbert space 
(corresponding to singlet dynamics) and black hole formation is expected to occur in the 
collapse of matter made out of non-singlet states. 

Note that we are not proposing that the symmetry of perturbative two dimen- 
sional string theory is broken non-perturbatively. Rather, we are saying that it does not 
prevent the appearance of non-trivial dynamics in the non-singlet sectors. As mentioned 
above, the situation is analogous to string theory on AdSs- There, the infinite dimensional 
spacetime Virasoro algebra (an analog of the symmetry of two dimensional string 
theory) is a symmetry of the full theory, but it does not prevent non-trivial dynamics and 
black hole formation in general, although the conformal block of the identity indeed does 
not contain black hole states (just like the singlet sector in two dimensional string theory). 

Finally, we would like to mention another aspect of our results, related to the Hage- 
dorn spectrum (|1.2j). In perturbative string theory, a spectrum of this sort would lead 



20 



For a review of gravitational physics in two dimensional string theory see [19|. 
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to an infrared instability (a tachyon), unless there was an almost complete cancellation 
between bosons and fermions [^6|j 47[|. It is an interesting open question whether this rela- 
tion between the density of states and IR stability is more general. Assuming that it is, it 
would be interesting to understand how to reconcile it with the behavior found here. 

One possibility is that the theory actually has fermionic excitations. Perturbative 
two dimensional string theory has no spacetime fermions; however, the statistics of the 
non-perturbative states is not obviously bosonic. Recall that in the singlet sector, the 
underlying degrees of freedom are free fermions. Another possibility is that the actual 
high energy density of states is in fact much smaller than that suggested by ( |1 . 2| ) . Further 
discussion of this possibility requires a better understanding of the high temperature phase 
of the theory (associated with tachyon condensation - see section 6). These issues deserve 
a better understanding. 
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Appendix A. Free compactified fermions from the MQM averaged over twist 
variable 



To show how free fermions emerge from the MQM partiton function ( |4.1|) let us use 
the following formula: 

Z S N ngl£t = J W] SU{ N)Z N {9) = JvM J [dn} SU (N)K 2nR (M,n+Mn), (A.l) 

where K2 n n(M, M') is the propagation kernel for the matrix variable from the initial value 
M to the final value M' during the time 2nR. 



Using the discrete time version of (|4.1| ) with the action given by 



S = J2 Tr (J^( M k- M k -i) 2 + aV(M k )^j , (A.2) 
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diagonalizing the matrices and applying consecutivly the Itzykson-Zuber formula to inte- 
grate out the angular variables |14[ we obtain in the limit of a small lattice spacing a — > 
the following representation for the propagation kernel of the eigenvalues u\,- ■ - ,itjv: 

K 2 -kr{{u}, {u'}) = 1 det kj K 27TR (u kl u'A (A.3) 
/\(u)A(u) J 

where A (it) = Ylk>j( u k ~ u j) ls the Vandermonde determinant. We ignored in ( |A.3| ) the 
overal normalization constant. Finally putting M = M' (which means u k = u' k ) and 
integrating with respect to the eigenvales of the final matrix M with the Dyson measure 
A 2 (u) Ylk=i we obtain from ( |A.1|) and ( |A.3| ) the following representation for the singlet 
partition function: 

r N 

Ringlet) = /JJ du . detkjK27TR ( Uk:Uj ), (A.4) 
^ i=l 

or, returning to the functional integral, 

N 

^singlet) = ^_ 1)V j TT Vu ^ e ~ f ^[U^f +V {u k )\ ^ 

where the sum goes over permutations V of the eigenvalues (u±, ■ • • , Wat) hi the process of 
propagation around the circle. 

The last formula gives the explicite free fermion representation of the singlet sector of 
the compactified MQM. 



Appendix B. Evaluation of the twisted partition function in the scaling limit 

If we forget for a moment about the wall, the problem decomposes into iV 2 independent 
inverted oscillators. They are related only by common twist angles. For each complex 
matrix element M^j, the Schrodinger equation on the wave function ty(Mkj) looks as 

"5 [mSmt, + M » Mt ) * (M *> ] = m{M ^ (B1) 

where there are no summations over the repited indices. We have to calculate from here 
the density of states in the presence of twists. In the lagrangean formulation we have the 
periodic boundary conditions Mkj(27rR) = e iekJ Mkj (0) , where 9kj = Ok — 9j. To express 
them in the hamiltonian language we have to go to the polar coordinates Mkj = rkj e l9kj , 
where we fixed the angular dependence by the difference of fixed twist angles, and look for 
the wave function in the form 

oo 

*{r kj ,d kj )= Xm(rkj)e im9k *. (B.2) 
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The density of energy will be given in terms of the sum of densities p m (E) of the states with 
a given angular momentum m weighted by e irndk: > = e im ( 9 k-0j) ( w fiich can be calculated 
from the analysis of the quasicalssical wave function, see |27j for the details): 

oo 

p(0,E)= eim9 Pm(E)- 

m= — oo 

The explicit calculation gives 

p( ( ,B) = !?f^f + i W logA». (B.3) 
smh[7T£/J sin 9 

Note that this result does not depend on the familiar cutoff of the c = 1 string theory for 
the nonequal twist angles 9 k j = 9 k — 9j 7^ 0. The system becomes stabilized due to the 
twisting! 

Now we can calculate the partition function of the individual matrix element in the 
fixed temperature T = 1/2-kR: 

r°° 1/2 
Z 1 (9 kj , R) = / dEe- 2 * RE p(6 kj , E) = > — . (B.4) 

J-oo COS(27T-R) - COS 9 kj 

Strictly speaking, this formula makes sense only for R < 1/2 and the distribution of 0's 
in the interval \9 max — d m in\ < kR- These limitations are inherited from the canonical 
(fixed N) ensemble used so far. We will see that it is not present in the grand canonical 
ensemble. 

From this formula, by taking the product of partition functions of the individual 
matrix elements, we derive the partition function of the twisted inverted oscillator 

Z N (e 1: ■■■,9 N ,R) = (2 sinnR)- N J[ 1 — (B.5) 

- L - L cos(27rit) — cos(6>i — 6j) 

where the first factor is the 9- independent contribution of the diagonal matrix elements. 



Appendix C. r-function and Hirota equations for the Toda chain hierarchy 

C.l. Vertex operator representation 

The simplest way (at least for a physicist) to obtain the fermionic representation is 
through the vertex operator construction. It is easy to see that this is the Coulomb gas 
partition function describing a chiral bosonic field with exponential interaction. Introduce 
the bosonic field (p(z) with mode expansion 

<p(z) =q + plogz + J2 —z~ n , (C.l) 

n=lO 
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[H n , H m ] = n5 m+n:0 ; [p, q] = 1. (C.2) 
and the vacuum state \l) denned by 

H n \l) = 0, (n>0); p|Z) = l\l). (C.3) 



The associated normal ordering is defined by putting H n ,n > to the right. Now define 
the vertex operators 

V q (z) =: e^ ql/2z) :: e -*(9~ 1/3 *) : (C.4) 

satisfying the OPE 



Then the r-function Q4.24p can be written as the expectation value 

ri[t] = (l | exp f J2 l nHn ) g exp I - f nHn ) | l) (C.6) 

\n>0 / \ n<0 / 

where the operator g is defined as 

g = exp (e 27rRfl I —V q (z)] . (C.7) 
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C.2. Fermionic representation 



The fermionic representation of the r-function now follows from the bosonization 
formulas 

if>(z) =: e~^ z) :, ip*(z) =: e^ z) :, d(p(z) =: ^*(z)V>(*) = (C.8) 

where 

V>(*) = £ W^^, V*(*)= E ^V - * (C9) 

is a chiral complex Ramond fermion field, whose modes satisfy the canonical anticommu- 
tation relations 

W> r ,tf]+ = 6 ra . (CIO) 



The Hamiltonians H n and the vertex operator ( |C.4|) are represented as fermion bilin 



ears 

H n = ■■^r-n^r: (^2) (C.ll) 

V q (z) = ^{q^ 2 z)t/j*(q- 1/2 z). (C.12) 
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and the vacuum states with given electric charge / satisfy 



<Z|VW = m* r = (r > /) 

VvlO=C r IO = o (r>l). 

C.3. Hirota bilinear equations 

The ensemble of r-functions with different charge satisfy a set of bilinear equations 
||48|| . Below we give a sketch of their derivation. 

The operator g defined by ( JC7| ) and ( |C12| ) is an exponential of a fermionic bilinear 
and as such can be thought of as a Gl(oo) rotation. The corresponding Lie algebra gl(oc) 
is generated by all bilinears J2 rs a r script ■ The operator 



S12 = V] A ® ip* = / ir~M 

^— ' 1 J 2m 



z) ®%l)*{z) 



is a fermionic analog of the tensor Casimir in the sense that it commutes with the tensor 
product of two copies of the operator g: 

S12 g®g = g®g 5i 2 . (C.14) 

Further, the definition of the vacuum implies that %j) r \ I) V) = (/ > /'), since either 

= or V) = 0. Therefore 

S 12 \ l)®\l') = 

and similarly for the left vacuum. 

Now we multiply eq. ( gHg) from the left by (l+l |e^»>o tnHn <g> (Z'-l | e E„ >0 t ™ w ™ an d 
..w,,, ,,,, > ™<o " n \l)®e 2^, n< o n ™ | /') and commute the fernion operators 

with the exponents until they hit the left (right) vacuum and disappear. As a result we 
obtain the following identities known as Hirota bilinear equations 



f %k*~ v exp ( ^ {tn - - ?+) Tt'V + c+) = 

"'C' 00 n>0 

/ ^"'' eX P(E(*-- £ nK n )^ + l(*-C-) ^_!(t' + C-). 

where 

C+ = (..., 0, 0, z" 1 , z~ 2 /2, z~ 3 /3, . . .), (- = (•••, * 3 A z 2 /2, z, 0, . . .). 



(C.15) 
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Expanding in y n = t' n — t n we obtain an infinite hierarchy of soliton partial differential 
equations. 



^2 P J+m (-2y+)Pj(D+) exp (^yjDjjn+m+ilt] ■ n[t] 

j=0 j^O 
oo 

^P i _ ro (-2y_)P i ( J D_)exp {^ Vj D>jr l+m [t] ■ r l+1 [t] 



(C.16) 



3=0 j^O 



where 



t± = (t±i, t± 2 , ...), y± = (y±i, y±a, •••)> 
£>± = (£>±i, £±2/2,^3/3,-) 
represent the Hirots's bilinear operators 



Dnf[t] ■ g[t] = ^-f(t n + x)g(t n - x) , (C.17) 

OX x=0 

and finally Pj are the Schur polynomials defined by 

00 00 

£>i[t] X J =exp(^t n A n ). (C.18) 

j=l n=l 

The lowest, Toda equation, which is what we will actually need, is obtained by equating 
to zero the coefficient in front of y-\. 

ndtd-m - dmd-iTt + n+M-i = 0. (C.19) 

Appendix D. Calculation of partition functions on torus and sphere 

D.l. Solution of the equation for the sphere 

Changing the variable y = e u and the function Xq = —(f) — 2u we bring the equation 
( fjHP to the form 

atifo = (d 2 u <j> + (d u <f>) 2 + 3d u <f> + 2)e^ 

where it does not contain explicitely u. 

By denoting ip{4>) = d u (p we obtain the first order ODE 

(ae-*- 1)^ = ^ + 1)^ + 2) 

with the solution 

(gu0±l) _ r 

where C is the integration constant. Solving the quadratic equation for d u (p and integrating 
it we obtain: 



4> + 2u = [ d(j> - 1 _ = - 2 arcsinh [ J 1 * Cc V */ 2 I + C 3 . 

J Wl-4Ca-4Ce* y/l - 4Ca \\ 4C J 

Going back to original variable y we see that fixing the couplings C = (^7^) 2 and C3 = 
we reproduce ( |5.7|) generating the series in Moore's formula ( |3.10| ). 
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D.2. Solution of the equation for the torus 

The easiest is to change the variable from y to t = e~ ~ x ° 1 according to the spherical 
solution (|5.7|). After a tedious but direct calculation ( |5.13| ) takes the form 

u ( t> i^i a f t ■ R ~ 1 h i n\ f i _ R ( R ~ 2 ) f l/(2-R) Cg - l ) 2t ~ 1 
[t-(R- l)]d t U + — (1 - l/t)f x - l2 t [{R _ 1)t _ ir 

It is a frst order inhomogenious ODE with the explicite solution 

R(R—2) f* R(s)-B(t) A/(2-m (i? - 1) 2 S - 1 



f , = R{R^- 2) J ^ eB{s) _ m s i/(2-i?) 



[{R-l)s-l]*[s-{R-l)] 



with 

B w =M f dt w^w\ = log ( <1/< " 2 ' [i " (fl " 1)1 

Fortunately (and may be due to the itegrability of the original Toda equation) the resulting 
integral for fi can be explicitly calculated: 

_ R(R-2) tV( 2 ~ fl ) / 2-R 1 1 

/l_ 12 t-(R-l) \2(R-1) [(R-l)t-l] 2 (R-l)t-l + 

where D is an integration constant. 

Using fi(y) = J dyf[ and integrating once more we get for the partition function 

, , , If 2DR R 2 \ , 

fM = -M{R^l + -(R^w) l ° gt 

1 f2DR 2 (R-2) 2R 2 -2R+1\, r/n . . 1 r „ ,„ ^, 4 

(D4) 

where the last term is the zero mode of Toda equation with yet unknown coefficient. 
To fix D and C let us use the boundary condition (|4.33| ) for A = 

/lOM) = - R+ 2 ^ log ft, A = 0. 

(Note that we always drop a constant, cutoff dependent, term from the partition function.) 
Since y — ► oo we have i = e^ 1-2 /^* — > y^ -2 and comparing the solution with the boudary 
conditions Q4.33| ) we conclude: 

2D ^^r + (R-i) 2 ~ R+R ■ 

Note that, "miraculously", the second logarithmic term disappears from ( p.l|) for this 
value of D. 

Since F(/j,, A) is regular at A = the constant C should be chosen as 

c _ R + R ~' 
12(2 — R)' 

With these values of C and D we finally obtain the result ( |5.14[ ). 
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